arXiv:l504.03525v3 [math.AP] 22 Mar 2016 


THE VARIABLE COEFFICIENT THIN OBSTACLE PROBLEM: 
OPTIMAL REGULARITY AND REGULARITY OF THE 
REGULAR FREE BOUNDARY 

HERBERT KOCH, ANGKANA RULAND, AND WENHUI SHI 


Abstract. This article deals with the variable coefficient thin obstacle prob¬ 
lem in n + 1 dimensions. We address the regular free boundary regularity, the 
behavior of the solution close to the free boundary and the optimal regularity 
of the solution in a low regularity set-up. 

We first discuss the case of zero obstacle and W 1,p metrics with pG (n +1, oo]. 
In this framework, we prove the C 1,a regularity of the regular free boundary 
and derive the leading order asymptotic expansion of solutions at regular free 

boundary points. We further show the optimal £' 1 > mm {i P 5 2 ^ regularity 
of solutions. New ingredients include the use of the Reifenberg flatness of the 
regular free boundary, the construction of an (almost) optimal barrier function 
and the introduction of an appropriate splitting of the solution. Important in¬ 
sights depend on the consideration of various intrinsic geometric structures. 
Based on variations of the arguments in [KRS15] and the present article, we 
then also discuss the case of non-zero and interior thin obstacles. We obtain 
the optimal regularity of the solutions and the regularity of the regular free 
boundary for W 1,p metrics and W 2,p obstacles with p E (2(n + 1), 00 ]. 
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1. Introduction 

In this article we continue our discussion of the variable coefficient thin obstacle 
or Signorini problem in a low regularity framework. Here our main objectives are 
an improved understanding of the (regular) free boundary, the determination of 
the asymptotic behavior of solutions close to the (regular) free boundary and the 
derivation of optimal regularity estimates for solutions. To achieve this in our low 
regularity set-up, we in particular introduce two key new arguments: The iden¬ 
tification of the regular free boundary as a Reifenberg flat set, which enables us 
to construct (almost) optimally scaling barrier functions (c.f. Section 3.1) and a 
“ splitting technique'" (c.f. Proposition 3.10), which allows us to deal with diver¬ 
gence form right hand sides and to identify the leading order contributions in the 
respective equations. 

Let us explain the precise set-up of our problem: We consider local minimizers 
of the constrained Dirichlet energy: 

( 1 ) J{ w ) = J (diw)(djiv)dx, 

Bf 

where we use the Einstein summation convention and assume that 

w G K ■.= {u G H\B+)\ v > 0 on B[ := B+ n {a; n+ i = 0}}. 

Here, the metric a®- 7 ' : B+ —> is a symmetric, uniformly elliptic tensor 

field which is W 1,p , p £ (n + 1, oo], regular and Bf := {x £ B\ C R" +1 | x n +\ > 0} 
denotes the upper half-ball. On the upper half-ball it is possible to consider arbi¬ 
trary local variations of local minimizers. Hence, local minimizers solve an elliptic 
divergence form equation in the upper half-ball, on this set they are “free”. How¬ 
ever, on the codimension one surface B[ they obey the convex constraint w > 0 
which leads to complementary or Signorini boundary conditions. In this sense the 
obstacle is “thin”. 

In the sequel we are in particular interested in obtaining an improved under¬ 
standing of the free boundary 

:= ds'^x e B [| w(x) > 0}. 
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This set (which for Lipschitz metrics a is of Hausdorff dimension n — 1, c.f. Remark 
3.1) separates the contact set , A„, := {x £ B[\ w(x) = 0}, in which the solution 
coincides with the obstacle, from the positivity set , Cl w := {x £ B[\ w(x) > 0}, in 
which the solution is “free”. Moreover, we seek to understand the structure of the 
solution close to the free boundary. 


Considering variations of local minimizers in the energy functional (1) leads to 
an equivalent formulation of the local minimization problem (1) in the form of a 
variational inequality posed in the energy space /C [ Jra87]: For a solution, w £ /C, 
of (1) we have 



v — w)dx > 0 for all v £ 1C. 


If in addition w £ H 2 (B+ ), this corresponds to an elliptic equation with comple¬ 
mentary or Signorini boundary conditions: 

dia^djW = 0 in Bt , 

( 2 ) 

w > 0, —a n+lj djW > 0, w(a n+1 ^djw) = 0 on B[. 

Here the Signorini condition is derived from the pointwise inequality —a n+1 d djiv(v— 
w) > 0 on B[ which holds for any v £ 1C. 


In the sequel we investigate the thin obstacle problem by studying solutions of 
(2). Moreover, we address variants of it which involve inhomogeneities, non-flat 
obstacles and boundaries and non-flat interior obstacles. 


1.1. Main results. In this article we first derive the C l,a regularity of the so-called 
regular set of the free boundary in the presence of W l,p metrics, a l fl with p > n+ 1. 
Moreover, in this framework we deduce a leading order aymptotic expansion of so¬ 
lutions to (2) with error estimates. Combining the regularity of the regular free 
boundary with the Carleman estimate from [KRS15], we then show the optimal 
£,i,mm{i —,?} re g U i ar ity of solutions with W 1,p metrics, for p > n+ 1. In 
addition to this, we also treat perturbations of the thin obstacle problem including 
non-flat free boundaries and obstacles, as well as inhomogeneities in the equations 
and the interior thin obstacle problem. 

In order to deduce these results, we rely on two main new ingredients: A “splitting 
argument” and the construction of (almost) optimally scaling barrier functions. 
The latter builds on the identification of the regular free boundary as a Reifenberg 
flat set. 

In the following subsections we elaborate on these results, put them into the context 
of the literature on the thin obstacle problem, and explain the main difficulties and 
the new arguments which are used to overcome these. 


We first recall the main results from [KRS15] on free boundary points: All free 
boundary points x £ T^, are classified by their associated vanishing order k x (c.f. 
Section 4 in [KRS15]): 

r™ = r 3 / 2 (wj) u (J r K (w). 

K> 2 
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Here r 3 / 2 (ro) is the so-called regular free boundary which is given by all free bound¬ 
ary points with vanishing order 3/2. The remaining set (J r K (u>) consists of all free 

k>2 

boundary points with a higher order of vanishing (c.f. Definition 4.1 in [KRS15] or 
Section 2.2 for the precise definition of k x ). Moreover, the map r„, E3 x K > k x is up¬ 
per semi-continuous. In [KRS15] we also proved that for each x £T W with k x < oo, 
there exists an L 2 -normalized blow-up sequence w x>r . such that the limit ui Xt o is 
a homogeneous global solution with homogeneity n x . Furthermore, if n x = |, the 
blow-up limit w x> o is two-dimensional and (up to a rotation of coordinates) is equal 
to c n Re(x n + ix n+1 ) 3/2 . 

In the first part of the present paper, we study the regular free boundary r 3 / 2 (u;), 
which is relatively open by the upper semi-continuity of n x . Relying on comparison 
principles (c.f. Proposition 3.17) combined with a splitting technique for equa¬ 
tions with divergence right hand sides (c.f. Proposition 3.10), we obtain the C x ' a 
regularity of the regular free boundary r 3 / 2 (w;): 

Theorem 1. Let a Riym 1 ' >x( ' n+1 ' ) be a uniformly elliptic, symmetric W 1,p , 
p £ (n + 1, oo], tensor field. Assume that w is a solution of the variable coefficient 
thin obstacle problem (2). For each xq £ r 3 / 2 (w), there exist a parameter a £ (0,1], 
a radius p = p(xo,w) and a C 1,a function g : B"(x o) —> R such that (possibly after 
a rotation) 

IR n B' p (x 0 ) = r 3 / 2 M n b' p {x 0 ) = {x| x n = g{x")} n b' p (x 0 ). 

We remark that we turn the usual order of the arguments around: Instead of 
first proving optimal regularity of the solution, and then regularity of the (regular) 
free boundary (c.f. [AC06], [ACS08], [PSU12]), we first prove regularity of the 
(regular) free boundary, and then deduce optimal regularity of the solution. The 
regular free boundary is C l a under conditions which do not imply C 1 ’? regularity 
of the solution. 

With the C 1,a regularity of the (regular) free boundary at hand, it then be¬ 
comes possible to study the local behavior of solutions around regular free bound¬ 
ary points. This is based on identifying the leading order in the asymptotics of 
solutions of (2) at regular free boundary points in the presence of W 1,p metrics 
with p £ {n + l,oo] (c.f. Proposition 4.6). As the asymptotics are complemented 
by a higher order error estimate, this allows us to obtain local growth bounds. 
Then combining this with our Carleman estimate, we are able to obtain the C ljl / 2 
optimal regularity of solutions associated with W 1,p metrics for p £ (2(n + l),oo]: 

Theorem 2 (Optimal regularity). Let aR : be a uniformly 

elliptic, symmetric W 1,p , p £ [2(n + l),oo], tensor field. Assume that ui is a 
solution of the variable coefficient thin obstacle problem (2). Then, there exists a 
constant C > 0 depending only on \\a 1 ^\\ w i, P ( B +y n,p such that 

IMI c i - i /2(b+ /2 ) < c IMI L 2( B +) • 

Finally, in the last part of the paper we prove that these results are not restricted 
to the flat thin obstacle problem, i.e. the setting in which the obstacle is the zero 
function. We show that, on the contrary, it is possible to deal with inhomogeneities , 
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non-constant obstacles, non-flat boundaries and even non-flat interior obstacles 
(c.f. Section 5). For instance, we prove the following result for non-flat obstacles: 

Theorem 3. Let a V : B+ —t be a uniformly elliptic, symmetric W 1,p , 

p £ (2(n+ l),oo], tensor field. Suppose that ip £ W 2,P (B' 1 ). Let w : B+ —> R be a 
solution of the thin obstacle problem 

difl 1 -’ dnW = 0 in Bt , 

(3) . 

w > tp, —a n+1 ^djW >0, (w — tp)(a n+1 ’i djiv) = 0 on B[. 

Then, the following statements hold: 

(i) The function w has the optimal Holder regularity: 

w e C^/ 2 {b+ /2 ). 

(ii) Assuming that 0 £ T^/ 2 (w), there exist a radius p > 0, a parameter a £ 
(0,1] and a C 1,a function g such that (potentially after a rotation) 

n B' p = r 3 / 2 (w) nfi'= {x\ x n = g(x")} n b' p . 

1.2. Literature and context. In the last years the thin obstacle problem and its 
variants have been a very active field of research. After the complete characteri¬ 
zation of the two-dimensional, constant coefficient thin obstacle problem by Lewy 
[Lew72] and Richardson [Ric78] as well as impressive partial results on the gen¬ 
eral problem [Fre77], [Kin81], [Ura85], a major new idea emerged only relatively 
recently: In [AC06] and [ACS08] Athanasopoulos, Caffarelli and Salsa introduced 
Almgren’s frequency function as a powerful tool of obtaining optimal regularity 
estimates and the C 1,a regularity of the free boundary for the constant coefficient 
problem. Here a is some constant in (0,1]. Later this was extended by [Sil07] and 
[CSS08] to the related obstacle problem for the fractional Laplacian. Moreover, in a 
recent article [KPS14] Koch, Petrosyan and Shi prove the analyticity of the regular 
free boundary for the constant coefficient operator by introducing a connection be¬ 
tween the thin obstacle problem and the Grushin Laplacian. Simultaneously, Savin 
and De Silva [DSS14] obtained the C°° regularity of the regular free boundary 
by exploiting higher order boundary Harnack inequalities. This discussion of the 
regularity of the regular free boundary is complemented by an article of Garofalo 
and Petrosyan [GP09] in which a monotonicity formula is used to characterize the 
structure of the singular set of the thin obstacle problem. 

While these results illustrate that there has been great progress in the constant 
coefficient thin obstacle problem, less is known in the variable coefficient framework. 
Here the work of Uraltseva [Ura85] has shown that it is possible to obtain C 1,a 
Holder regularity of solutions of the thin obstacle problem in the presence of bF 1,p , 
p £ (n + 1, oo], metrics. In her result the Holder exponent is some value a £ (0,1/2] 
which depends on the ellipticity constants of the coefficients and the value of p. 
Only very recently, the variable coefficient problem has been further investigated: 
In [Gui09] Guillen deals with the problem in the context of Holder regular C 1,a , for 
some a £ (0,1), metrics. This is improved in an article by Garofalo and Smit Vega 
Garcia [GSVG14] in which the authors derive the optimal regularity of solutions 
of the thin obstacle problem in the presence of C 0,1 metrics and C 1 ’ 1 obstacles. 
This argument is based on an extension of Almgren’s monotonicity formula to 
the setting of low regularity metrics and obstacles. In a recent preprint [GPG15] 
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Garofalo, Petrosyan and Smit Vega Garcia further build on this and deduce the 
C 1,a , for some a £ (0,1), regularity of the regular free boundary in the presence of 
C 0 ' 1 metrics and C 1,1 obstacles by proving an epiperimetric inequality. 
Complementing these results by relying on Carleman estimates instead of frequency 
formulae and working in the setting of Sobolev metrics a lJ £ W 1,p , p £ (n + 1, oo], 
[KRS15] provides an alternative proof of the almost optimal regularity of solutions 
to (2). In the present article we extend these ideas, and, building on our previous 
work, prove optimal regularity in the framework of W 1,p , p £ (n + l,oo], metrics 
as well as the C 1,a regularity of the regular free boundary. 

1.3. Difficulties and strategy. As already in [KRS15] the central difficulty with 
which we deal throughout the article is the low regularity of the metric. 

Building on the results from [KRS15], we analyze the (regular) free boundary 
(c.f. Section 3). Traditionally, it is studied by relying on comparison principles 
(c.f. [PSU12]): Differentiating the equation (2), an analysis of the equation for the 
tangential derivatives allows to transfer positivity properties of derivatives of the 
blow-up solutions to the problem at hand. This then permits to conclude the C 1,a 
regularity of the regular free boundary. 

In a low regularity set-up this becomes more difficult. In particular we have to 
include divergence form right hand sides in our discussion. Hence, as a key new 
ingredient we exploit a “splitting technique” in which we divide our solution into a 
“controlled error” which handles the low regularity contributions originating from 
the metric, and a “main part” which captures the behavior of our solutions (c.f. 
Proposition 3.10). This allows us to argue along similar lines as in the literature. 
Yet, we have to overcome a second difficulty: In order to provide a framework which 
also deals with non-flat obstacles of low regularity, we have to work as close as pos¬ 
sible to the scaling critical setting. Thus, instead of proving linear non-degeneracy 
of the tangential derivatives, we prove a (nearly) square root non-degeneracy in 
appropriate cones (c.f. Proposition 3.17). In order to achieve this, we introduce 
a second main new ingredient and construct a new barrier function which exploits 
the Reifenberg flatness of the free boundary (c.f. Proposition 3.16). 

In the second part of our argument, we return to the investigation of solutions 
of the thin obstacle problem (c.f. Section 4). Here we rely on the free boundary 
regularity which allows to improve the almost optimal growth estimates which were 
previously obtained from the Carleman estimate (c.f. Lemma 4.1 in [KRS15] and 
Corollary 4.11). 

We argue in two steps in which we combine local with global information: First 
we prove a local growth estimate around regular free boundary points in which we 
obtain the optimal growth estimates. As these bounds however rely on comparison 
arguments which depend on the free boundary itself, they are not uniform in the 
free boundary points. By virtue of the C l,a regularity of the (regular) free bound¬ 
ary, these growth estimates can be obtained by an asymptotic expansion of solutions 
around regular free boundary points. For the identification of the leading order con¬ 
tribution of the expansion and of the corresponding higher order error estimates, 
we exploit our boundary Harnack estimate in combination with the boundary reg¬ 
ularity for W 1,p metrics with p £ (n+ 1, oo]. Here we exploit our main new results 
on the free boundary regularity. 
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In the second step in our optimal regularity argument for W 1,p metrics with p > 
2 {n + 1), we then combine the optimal, but non-uniform local with optimal global 
information. For this we rely on Lemma 4.2 in [KRS15] which allows us to transfer 
the local into global information. This is the only point at which we directly return 
to the arguments from [KRS15]. 

Finally, in Section 5 we comment on the stability of the described methods by 
applying them to variants of the thin obstacle problem. Using the scaling of the 
Carleman inequality, we first show that it is possible to deal with inhomogeneities 
(c.f. Section 5.1). This then immediately entails that all the previous results remain 
true for sufficiently regular, non-constant obstacles, although we only require that 
they are W 2,p regular for some p > 2(n + 1) (c.f. Section 5.2). 

Concluding the section on variants of the thin obstacle problem, we discuss the 
setting of the interior thin obstacle problem (c.f. Section 5.4). Here we are con¬ 
fronted with additional difficulties, which arise due to a slightly modified boundary 
condition (instead of a sign condition on the Neumann derivative, there is a sign 
condition on the fluxes across the interior boundary). This leads to slight mod¬ 
ifications in the derivation of the Carleman inequality from [KRS15] and yields 
leading order linear contributions in the asymptotics of the Neumann derivative 
(c.f. Proposition 5.6). 

1.4. Organization of the paper. After briefly recalling auxiliary results, conven¬ 
tions and the notation from [KRS15] in Section 2.2, we begin with the analysis of 
the regular free boundary in Section 3. Here Proposition 3.10 is a crucial technical 
tool to overcome the difficulties with the low regularity set up. In Section 3.2 we 
first prove the Lipschitz regularity of the regular free boundary (c.f. Proposition 
3.8). This is based on the observation that the free boundary is Reifenberg flat, 
which we exploit in Section 3.2.1, in order to construct an appropriate, sufficiently 
well scaling barrier function (c.f. Proposition 3.16). In Section 3.3 we improve the 
Lipschitz regularity to gain C 1,a regularity (c.f. Proposition 3.22). 

In the second part of the paper we return to the study of solutions of (4): In Sec¬ 
tion 4 we identify the leading order asymptotics of solutions of the thin obstacle 
problem (Proposition 4.3). This allows to derive growth bounds (Corollary 4.9 and 
Corollary 4.11) as well as the optimal regularity of solutions of (2) (c.f. Theorem 4). 
Finally, in the last part of the article, in Section 5, we illustrate how the previous 
results can be transferred to variations of the thin obstacle problem. 

2. Preliminaries 

In this section, we briefly explain our normalizations and notational conventions. 

2.1. Auxiliary results. We start by recalling that, due to the discussion in Section 
2.1 in [KRS15], we may, without loss of generality, consider solutions w of (2) with 

(AO) IMIl2(b+(o)) = 1, 

(Al) a l,n+1 {x ', 0) = 0 on M" x {0} for i = 1,..., n, 

(A2) a 10 is symmetric and uniformly elliptic with eigenvalues in the interval 

[ 1 / 2 , 2 ]. 

In addition, we in the sequel also make the following hypotheses (where we either 
assume (A3) or (A3’)): 
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(A3) a IJ £ W 1 ’ P (B i (0)) for some p £ (n + 1, oo], 

(A3’) a lJ £ W 1,p (Bi (0)) for some p £ (2 (n + 1), oo], 

(A4) a ij {0)=6 ij . 

These assumptions allow us to reduce (2) to 

dia^dnW = 0 in Bf , 

(4) 3 \ 

w > 0,—dn+iw > 0, w(d n +iw) = 0 on B[. 

Due to the H 2 estimates and the almost optimal regularity result of [KRS15], it 
is possible to interpret (4) not only as a variational inequality in the energy space 
H 1 (B j*~), but also to understand the equation and its boundary values in a classical 
pointwise sense. 

Let us comment on the conditions (A0)-(A4): We recall that it is always possible 
to achieve (AO) by a suitable normalization. Condition (Al) is a consequence of an 
appropriate change of coordinates, c.f. Uraltseva [Ura85]. Also, assumption (A2) 
and (A4) can always be achieved by an additional affine change of coordinates (and 
a rescaling) and thus do not pose additional restrictions on our set-up. 

Conditions (A3) and (A3’) are regularity assumptions which allow us to apply the 
results of [KRS15]. Here condition (A3’) is slightly more restrictive. While this 
assumption is not needed in the case of flat obstacles, it becomes necessary in our 
treatment of non-flat obstacles, as a consequence of our strategy of proof: We work 
on the level of the differentiated equation. Moreover, we stress that the integrability 
assumption p > 2 (n + 1) yields the embedding W 1,p C 0,1 / 2 which, by interior 
regularity, is needed (and might also suffice) to derive the (optimal) C 1,1 / 2 regularity 
of solutions to the variable coefficient thin obstacle problem. Both conditions (A3) 
and (A3’) imply that 

| a ij (x) - 5 lj \ < C n ||Va ij || iI , (jB , ) for all x £ B+ 

by Morrey’s inequality. 

2.2. Notation. We use the same notation as in [KRS15], which we briefly recall 
in the sequel: 

• R” +1 := {x £ R n+1 | x n+1 > 0}, IT +1 := {x £ R n+1 | x n+1 < 0}. 

• For points x £ R ra+1 we also use the notation x = (a/, x„+i) or x = 
(x", x n , x n +i), if we want to emphasize the roles of the respective lower 
dimensional coordinates. 

• Let xo = (xq,0) £ R” +1 . For the upper half-ball of radius r > 0 around 
Xq we write Bfl(x o) := {x £ R” + | \x — Xo| < r jS the projection onto the 
boundary of R" +1 is respectively denoted by B' r (x o) := {x £ M n x {0}| |x — 
Xo| < r}. If xo = (0,0) we also write Bfl and B' r . Analogous conventions 
are used for balls in the lower half sphere: B~(x o). Moreover, we use the 
notation B"(xq) = {x £ M n_1 x {(xo) n } x {0}| \x" — x'q\ < r}, where 
xo = (xq, (x o )n,0),x = (x",x n ,0). 

• Annuli around a point Xq = (xq,0) in the upper half-space with radii 0 < 
r < R < oo as well as their projections onto the boundary of R" +1 are 
denoted by by A+ fl (x 0 ) := B^(x 0 ) \ B+(x 0 ) and A' r R (x 0 ) := B' f { (x 0 ) \ 
B' r (x o) respectively. For annuli around xo = (0,0) we also omit the center 
point. Furthermore, we set A t ,r{x o) := A+ (x 0 ) U A r R (x 0 ). 
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We use C v (e n ) to denote an (n + l)-dimensional cone with with opening 
angle rj and axis e„. Analogously, C^(e„) refers to a flat (i.e. n-dimensional) 
cone on {x n+ i = 0} with opening angle 77 and axis e n . 

For / : R n+1 -4 Rwe set V'/ and V"/ to denote the derivatives with 
respect to the x' and x" components of x. 

We use | • | to denote the standard Euclidean norm. 

Distances with respect to a Riemannian metric g (e.g. if they are induced 
by certain operators as in Section 3.4) are denoted by d g (-, •). 
dist_y(X, Y) := maxjsup d(x, Y), sup d(y, X)} denotes the Hausdorff dis- 

x£X y£Y 


tance of two subsets X and Y in W n . 

Let w : B+ —> 1 be a solution of (4). Then 

— Q w := {1 e I" x {0} | w(x) > 0} denotes the positivity set. 

— := ds'^w is the free boundary. 

— A w := B[ \ f l w is the coincidence set. 

— T3{w) := {x € T w \k x = |} C is the regular set or the regular free 

lnCr- 2 ^ 1 IMIj-a^-f. {x)) ) 

boundary. Here k x := limsup-r—pj—- is the vanishing 


order of u> at ax 

As in [ [RS15] we use the notation w r (x) := „ + i u, ( rx -- to denote 

the L 2 normalized rescaling of w at zero. Analogously we consider the L 2 - 
normalized blow-ups at arbitrary free boundary points xq £ r^: w X0:r (x) := 

w(xo~\-rx) 

r-^\M\ L 2 (B + (xo)) 

We set w 3 / 2 (x) := c n Re(x„ + ix n+ 1 ) 3 / 2 , where c n is chosen such that 


\W3/2\ 


= 1. 


H 2 (B+) 

In the following we often use the abbreviation := (2 y/n) -1 as well as 

C * := ILp(Bi) - 

We use Lq = a^dij and L = dia^dj to denote the non-divergence and 
divergence form operators involved. 

We reserve the parameter e > 0 to quantify the Reifenberg flatness of a 
domain T C R" x {0}, as well as the parameter eo to measure the closeness 
of a solution of (4) to the L 2 normalized model solution w 3 / 2 {x). 

We use the notation A < B to denote that there exists an only dimension 
dependent constant such that A < CB. Similar conventions are used for 


> 


3. Regularity of the Regular Free Boundary 

In this section we deduce the C 1,a regularity of the regular part of the free 
boundary for solutions of the thin obstacle problem (4) and thus provide the proof 
of Theorem 1. 

We recall that it is possible to classify free boundary points by their vanishing 
order k (Section 4 of [KRS15]): 

r«, = r|(w) u (J r K M, 

k>2 
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where V K (w) := {x £ T™] n x = k}. The set T|(w;) is referred to as the regular set 
of w or the regular free boundary of w. A corresponding point x £ T^w) is called 
regular. 

Although we are mainly interested in the regular free boundary, in passing, 
we make the following observation on the Hausdorff-dimension of the (whole) free 
boundary: 

Remark 3.1 (Hausdorff dimension of T^,). We claim that, if the obstacle is zero 
and the metric a lJ is Lipschitz continuous, the Hausdorff-dimension of the (whole) 
free boundary is less than or equal to n — 1. Indeed, this follows as in Lemma 
).l and Theorem 8.10 in [SW]. Similarly as in these results, we prove a slightly 
stronger statement: We show that 

H s (S w D B[) = 0 for all s > n — 1, 

where S w := {x £ B(f \ w(x) = |V?w(x)| = 0} in particular contains T^,. 

We only provide a sketch of the proof here and refer to [SW] and [Sim96], Chapter 
3.3 for the details. As in [SW] we argue in two steps and first show that ifip £ C 1,a 
for some a > 0 is a solution of the constant coefficient thin obstacle problem, then 
the Hausdorff dimension of the associated free boundary S v fl B[ is less than or 
equal to n — 1. 

The argument for this follows by contradiction: Assuming that the statement of the 
claim were wrong, there existed s > n — 1 such that 

TL s (S,p flB[) > 0. 

By density arguments for the Hausdorff measure (c.f. [EG15], Chapter 2.3, Theo¬ 
rem 2.7) there exists a point z £ S v fl B[ such that 

6*^(S v ,z) := limsupp~ s n s (S v n B' p (z)) > 0, 

p—>• 0 

oo oo 

where for A C B[ we have p s (A) := inf Y) P S j an d A C (J B' (yA and yj £ B[. 

3=1 3=1 

Hence there exists a sequence <jj > 0 with aj 0 such that 

lim vj s Hs(S v n B' (z)) > 0. 

<Tj—>0 J 3 

Rescaling the function by setting Pj(x) := _ - wji/j x £ ^ cr - 1 (i_| z |) 

results in 

liminf ti s {S Vj D B[) >0. 

By compactness (c.f. [SW], [ - > STJ12] y ) we have that along a subsequence 

<Pj> in C( oc (B+), 

where if £ C 1,a fl W 2,2 is a homogeneous solution with degree larger than or equal 
to 1 + a and it is normalized such that ||V , IIl 2 ( B +) = 1. Moreover, as in [SV ] 

n B[) > 0. 

We now repeat the outlined argument with if instead of tp. This yields a function 
ifi with p. s (S^, 1 nB]) >0, 0 £ and ip i being invariant under composition with 
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translations in the direction A z for all A £ 1. Indeed, the last point follows from 
the observation that J\f^ 1 (0) =J\f^ 1 (z), where 

P f \Difi{y)\ 2 dy 
■M/>! (x) := lim BpiX) — 

J m(y)[dy 

dB p (x) 

denotes the frequency function (c.f. [SW], Section 2 and [PSU12], Chapter 9 for 
the existence of this limit and its properties). We also note that this is well-defined 
for x £ and not only for x £ F^. Similarly as in Remark 2.) (2) in [SW] this 
implies that for the homogeneous homogeneous function ipi it holds ipi (x + A z) = 
ifi(x) for all x £ M” +1 . Arguing inductively and carrying out the argument a further 
n — 1 times (which is possible by our contradiction argument), we end up with a 
function if n (x i,... ,x n +i) = f(x n + 1 ) which only depends on the x n +\-variable and 
is non-trivial. As it is however also harmonic and satisfies if n ( 0) = |V^ n (0)| = 0, 
this yields a contradiction and hence results in our claim. 

In a second step, the more general case of a solution w to the variable coefficient 
problem (4) with Lipschitz coefficients a lJ is treated. Here the claim is that also in 
this situation the Hausdorff-dimension of S w D By 2 is less than or equal to n— 1, 
i.e. for all s > n — 1 it holds H. S (S W f~l By 2 ) = 0. Again this argument follows by 
contradiction assuming the claim on the Hausdorff dimension were wrong. Then 
however blow-up argument (in combination with compactness, which holds in the 
setting of Lipschitz coefficients, c.f. for instance [S\ }) reduces the situation to 
that of a non-trivial solution to the constant coefficient problem with Hausdorff 
dimension larger than n — 1. This however is in contradiction with our first claim. 

Returning to our main problem, in the following we study the regularity prop¬ 
erties of the regular free boundary F|(«;). This is divided into several steps: In 
Section 3.1 we begin by proving that the regular free boundary is a relatively open 
set of the free boundary and that, due to the good compactness properties of blow¬ 
up solutions, it is Reifenberg flat. Then in the following Section 3.2 we deduce 
its Lipschitz regularity. Our argument for this strongly relies on the boundary’s 
Reifenberg flatness, as this allows to construct (nearly) optimally scaling lower bar¬ 
rier functions. In Section 3.3 we exploit the Lipschitz regularity to infer the C 1,a 
regularity of the regular free boundary. Here we argue via a Carleson and a bound¬ 
ary Harnack estimate. Last but not least, in Section 3.4 we provide the proof of 
one of our most central technical tools, the splitting technique which is stated in 
Proposition 3.10. 

3.1. Reifenberg flatness. We begin the investigation of the regular set by showing 
that it is a relatively open subset of : 

Proposition 3.2 (Relative openness of r 3 / 2 (w;)). Let : R+ —> X < ' n+l ^ be 

a uniformly elliptic W 1,p tensor field with p £ (n + l,oo]. Let w be a solution of 
(4) in B(~. Then r 3 / 2 (w;) is a relatively open subset ofT w . 

Proof. By Proposition 4.2 in [KRS15] the mapping T,,, 9 x 1-9 k x is upper semi- 
continuous. Therefore, the preimage of the set of free boundary points with n x < 2 
is relatively open. Due to the classification of the possible homogeneous solutions 
and due to the regularity of general solutions, k x < 2 already implies that n x = i 
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(c.f. Proposition 4.6 in [ RS15]). Hence, we infer that rj(u>) is a relatively open 
set of the free boundary. □ 

In order to study the regular set in greater detail, we recall a strong compactness 
result for L 2 normalized blow-ups of solutions of (4). In contrast to the setting of 
free boundary points of higher vanishing order, it is possible to show that at regular 
free boundary points any blow-up is a 3/2 -homogeneous global solution. 

For convenience of notation, in the sequel we set 

(5) w 3 / 2 {x) := c n Re(x n + ix n+ i) 3/2 for x G Bf , 

where c n > 0 is an only dimension dependent constant, which is chosen such that 

II u; 3 / 2 IL 2 (b+) = 

Using this convention, we formulate the L 2 normalized blow-up result at regular 
free boundary points: 

Lemma 3.3. Let abe a uniformly elliptic W 1,p tensor field, 
with p G (n +1, oo]. Let w be a solution of (4) in Bf. Consider xo g r 3/2 (w)n5; /2 

and w rXo {x) := _ n _ | _ 1 u ’U°+ r:r ) -. Then for any sequence with rj —> 0 

r__a_IHI L 2 lB + lmo)) 

there exists a subsequence {rj k }ken and a matrix Q G SO{n + 1) such that 

( 6 ) w rjk >Xo (x) -* w 3 / 2 (Qx) in C 1 (B^ 2 ). 

Proof. The proof of (6) relies on an argument by Andersson [Andl3] and reduces 
the problem to Benedicks’s theorem [Ben80]. For completeness we give a brief 
sketch of it. 

To this end, we first note that by Proposition 4.4 in [KRS15] there exists a subse¬ 
quence {Tj fe }fcgN such that 

w r jk ,x 0 -> wo in C 1 (B+ /2 ). 

Here wq is a global solution of the constant coefficient thin obstacle problem with 
II u 'oIIl 2 (b + ) = 1- By Remark 17 in [KRS15] we further infer that wq grows of a rate 
R K with k G (1,2) at infinity. Analogously, any tangential derivative d e wo , with 
e G S n C\ B[, has at most a growth rate of i? K_1 with k G (1,2) at infinity. 

Next we show that there exists Q G SO{n+ 1) such that wq{x) = w 3 / 2 {Qx). Indeed, 
due to the at most R K_1 , k G (1,2), growth at infinity of d e wo, the Friedland- 
Hayman inequality [FH76] (see also Lemma 5.1 in [Andl3]), yields that for any 
tangential vector e G S n (~l B[ the directional derivative d e wo has a fixed sign in 
R" +1 . Thus, arguing as in Proposition 9.9 in [PSU12], we have that w o only depends 
on two variables, which, up to a rotation Q, can be written as wo(x n ,x n +i). As a 
consequence of the sign condition on d n wo (without loss of generality we assume 
d n wo > 0), the fact that <9 n+ iw;(0) = 0 (due to the C 1 convergence), and the 
growth rate at infinity, the coincidence set Ao of wq is the n-dimensional half-plane 
Ao = {x n < Q,x n +i = 0}. Due to the Signorini condition, after an even reflection 
about x n +i, d n wo solves 


Ad n wo = 0 in K” +1 \ A 0 , d n wo = 0 on A 0 , d n wo > 0 in R” +1 \ A 0 . 
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Applying Benedicks’s theorem (c.f. Proposition 1 in [Andl3] and also [Ben80]) to 
d n w o (or a Liouville type theorem for harmonic functions after opening up the do¬ 
main), we have that up to a multiplicative constant d n wo{x) = Re(a; r j + *|a; IJ ,+i|) 1 / 2 . 
A similar argument applied to d n +iWo yields d n +iivo(x) = lm.(x n + ix n + i ) 1 / 2 in 
B+. Hence, necessarily vjq{x) = c n Re(a; n + ix n +i) 3 ^ 2 = 1 / 13 / 2 ( 2 ;) in B+ (where we 
used that ||ui 0 || i 2 (B + ) = IK/ 2 || L 2 (B +) = !)• □ 

Keeping this compactness result in the back of our minds, we proceed by showing 
that the regular free boundary is Reifenberg flat. For this we first recall the following 
definition (c.f. [CKL05], [LMS12]), which allows us to infer our first regularity 
result, Proposition 3.5, for the free boundary. 

Definition 3.4 (Reifenberg flatness). A locally compact set T C K™ is ( 6,R ) 
Reifenberg flat if for every xq £ T and every r £ (0,R] there is a hyperplane 
L(r,xo) containing Xq such that 

(7) - distir (r n B r (x 0 ), L(r, x 0 ) n B r {x 0 )) < 5 , 
r 

where dist/r(X, Y) := maxjsup d(x, Y), sup d(y, A)} denotes the Hausdorff dis- 

xeX y£Y 

tance of two subsets X and Y in K. m . 

Recalling the convention from (5), we have the following regularity result: 

Proposition 3.5. Let a u : B+ —> be a uniformly elliptic W 1,P {B^) 

tensor field with p £ (n + l,oo] which satisfies (Al), (A2) and (Af). Let w be a 
solution of (4) in B±. Then for each e > 0 there exists 5 > 0 such that if 

(i) ||ui- W3/2|| c?1(B + ) < S, 

(ii) II^-^IL=c(b+| <^ 

then for each Xq £ D By 2 and each radius r £ (0,1/2), there is a rotation 
S = S(r, Xo) £ SO(n + 1) such that 

(8) ||«V,x o (20 - w 3 / 2 (Sx)\\ c1{b +j < e. 

In particular, if S is chosen sufficiently small, then r^, fl B^ 2 = r 3 / 2 (w>) fl R(y 2 . 

Proof. We argue in two steps. First, we derive a non-degeneracy condition at free 
boundary points, which is crucial in order to obtain good compactness properties. 
This first step also immediately implies that r 3 / 2 (ui) fl B[/ 2 = F^, fl R(y 2 . Then, in 
the second step, we argue that ( 8 ) holds. 


Step 1: Compactness. We begin by proving the following claim: 


Claim. Assume that the conditions of Proposition 3.5 are satisfied. Then for each 
A > 0, if 5 = S(n,p, A) is chosen sufficiently small, for all r £ (0,1/2) and for all 
Xq £ r w H B' 1/2 , we have that 


(9) 


n + 1 

r 2 


IMIl 2 (B+(x 0 )) — 


1 a _i_\ 

— 7 * 2 

2 


We prove the claim as a consequence of the closeness assumption (i) to the model 
solution and of Corollary 3.1 in [KRS15] for a suitable choice of r. 
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Consider R := (4(5) n + 4 with 0 < <5 < (1) 2+4 . Then, by the triangle inequality 
and by condition (i), we have that for any xo G T,„ D B^ 2 

IMIl 2 ( B+(* o)) - ll W 3/2|| i 2 (B + (a;o)) ~S> \\w 3/2 \\ l2{b +^- o)) - S. 

for all r G [i?, 1). Here Xq G T,„ 3/2 D denotes a point such that \xo — xo\ < 6 
(which exists by our assumption (i)). By virtue of the scaling of W 3/2 at free 
boundary points and by the choice of R, we thus infer that for all r G [R, 1/2) 

( 10 ) r ~^ IMIz^b+^o)) ^ \ ri - 

This shows (9) for all r G [R, 1/2). 

We now argue that (9) is also true for r G (0,1?): Let r G (0,1?) and r = 

1+( 1 . |iL — ^ 2< hi(r) o)) ) -— (where cq is the constant in Corollary 3.1 in [KRS15]). 

Applying Corollary 3.1 in [KRS15] at the scales r, R and 1/2 (where R is defined 
as above) yields: 


e r0(lnR)??— l| ln ??|—2|| w || 


L 2 (A+ 2 R (^o)) ^ 


C(n,p) ( 


g r0( Inr) r ~l 


HI 


L 2 (B+(x 0)) 


+e r0(ln(l/4)) 


I L*(B+ /2 (x 0 )) J 

< C{n,p). 

Here the last line follows from the definition of r. As a result, 

II w IIl 2 (a+ 2R (x 0 )) ^ C(n,p)e^ lnR) R\lnR\ 2 . 

Combining this with (10) and that (j){t) = —(1 + cof )t — cq (l + ln(— t) + 0( j)) as 
t — 00 , we infer that for our choice of R = (45) "+4 

1 l ’ ( a +„(*„)) < C(n,p)R^+^eV^R\lnR\^. 

Hence, using that Cq, t < 1 


\r^ < Ml 


7 r 

1 + Co 2 


Ccr 


| hr l?f) 


T + 1 < 


ln(C(n,p)) 4 ln(|lnl?|) 


2 In (1?) In (1?) 

which, by plugging in the expression for r, can be rewritten as 
Ccq 


1 + 


> 


(1 + c 0 f)|lnl?| 
n + 4 In (C(n,p)) 


MIM| L 2 (b + (xo)) ) 

In | ln(l?)| 


-4- 


Ccn 


ln(r). 


2 In (1?) In (R) l + c 0 ^\lnR\ / 

Thus, a sufficiently small choice of R (and thus S depending on A) then yields (9). 
Recalling that r G (0, R) was chosen arbitrarily and combining the just derived 
estimate with (10), finally yields the full non-degeneracy condition (9) (i.e. for the 
full range of radii r G (0,1/2)) and hence concludes the proof of the claim. 

As a direct consequence of the claim, we note that any free boundary point 
Xq G r w n B[/ 2 has vanishing order at most § + A. By choosing A < we infer that 
k Xq G [0,2). By Proposition 4.6 in [KRS15] this implies that k Xo = 3/2. Therefore, 
T w nB' 1 / 2 =T 3 / 2 (w)nB' 1/2 . 
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Step 2: Proof of (8). We argue by contradiction. Suppose that (8) were wrong. 
Then there existed a parameter eo > 0, a sequence 8 k —> 0 and a sequence, w k , of 
solutions to the thin obstacle problem with associated metrics d ?, i.e. 

didj!djW k = 0 in B+, 
w k > 0, — d n+ \w k > 0, w k d n+ \w k = 0 on B [, 

as well as a sequence of radii r k £ (0,1/2) and a sequence of points Xk £ C\B[^ 2 

such that 

(i) \\w k - w 3/2 \\ c i {B +)< s k, 

(ii) bk -^ j || L =o( B +) < fa, 

but for all S £ SO(n + 1) 

(U) \\ W r k ,x k ( X ) ~ W 3/2( Sx )\\ C i( B + 2 ) ^ £ 0' 

We can also assume that fl B[j 2 9 Xk —i► Xo £ T W3/2 (~l B^ 2 and that —>■ r q. 

We now show how the compactness statement of step 1 yields a contradiction to 
our assumption (11). For this we consider the L 2 normalized rescalings, w k x (x), 
of w at a free boundary point x £ T^ie D B[/ 2 - As the conditions of the claim 
from step 1 are satisfied at any free boundary point x £ fl the doubling 
property (c.f. Proposition 4.3 in [KRS15]) is satisfied uniformly in the choice of 
the free boundary point x £ and in the sequence w k . Hence, considering the 
sequences rk, Xk from above, entails that on the one hand w k kX Ax) is bounded in 
H 1 (B^~). Therefore, up to a subsequence it converges strongly in L 2 (B+) to a limit 
wq(x) £ H 1 {B+) with ||wo|| L2 ( B +) = 1. By the L 2 — C 1,a estimates this also entails 
that w k k Xk (x) -£ vjq in C ,1 (i?)y 2 ). On the other hand, we recall the C 1 convergence 

(12) w k {x) -A w 3 / 2 (x), 

which follows from (i). We now distinguish two cases: 

Case 1: r 0 > 0. In this case we observe that due to (12) 

w k (x k + r k x) u> 3 / 2 (xo + r 0 x) = rg /2 w 3/2 (x), 

IhIL »(*+(*>» ^ r o“^ \M\ L2{<{xo)) =r 3 0 /2 , 

where in the last line we used the L 2 normalization of w 3 / 2 and that 
w 3 / 2 is homogeneous of degree 3/2. Consequently, w k Xk (x) —> w 3 / 2 (x) 
in C 1 )!!^), which is a contradiction to (11). 

Case 2: ro = 0. In this case we note that, since w k is a solution of the thin obstacle 
problem, the same is true for w k Ax). Moreover, by assumption (ii) 
ajf —> S lJ , which implies that wo(a:) is a nontrivial global solution of the 
constant coefficient thin obstacle problem which satisfies = 1. 

By (9) and the uniform upper growth estimate of Lemma 4.1 in [KRS15], 
the growth of wo at infinity is of the rate with n £ (1,2). Thus, by the 
proof of Lemma 3.3, w 3 (x) = w 3 / 2 (Sx) for some rotation S £ SO{n + 1). 
This yields a contradiction to (11). 

Combining the previous two cases shows (8). □ 
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Proposition 3.5 immediately entails the following Corollary: 

Corollary 3.6. Let a*-' : B /" —y be a uniformly elliptic W 1,P (B /“) 

tensor field with p £ (n + l,oo] which satisfies (Al), (A2) and (A4). Let w be a 
solution of (f) in B± . For each e > 0 there exists S > 0 such that if 

(i) \\w - w 3/2 \\ cl(B + ) < 6, 

(ii) ~(b+)<*> 

then the free boundary r w D * s (c 1/2) Reifenberg flat. 

Proof. We regard P,,, D B[^ 2 as a subset of K ra x {0}. By Proposition 3.5, we have 
that for each r £ (0,1/2) and Xo £ fl B^ 2 , there exists S(r,x 0 ) £ SO(n + 1) 
with 

(13) ||w a:o ,r(S , (r,xo)x) - w 3/2 {x)\\ cl{B + ) < e. 

Consider L(r,x o) := (xo + x £ K” x {0}| x ■ {S(r,x o)e n ) = 0}. Due to (13), we 
obtain that 

- dist#(r„, n B r (x o), L(r, x 0 ) D B r (x 0 )) < e, 
r 

which shows the (e, 1/2) Reifenberg flatness of Tu, fl B' x 2 . □ 

By using an L 2 normalized blow-up, this can be transferred to obtain the Reifen¬ 
berg flatness of the regular free boundary of a general solution of the thin obstacle 
problem (which does not necessarily satisfy the closeness assumption (i)): 

Corollary 3.7 (Reifenberg flatness of r3/ 2 (u>)). Let a y : B+ —>• Riym 1 ^ x ^" +1) be a 
uniformly elliptic W 1,p tensor field with p £ (n+l,oo]. Letw be a solution of (4) in 
B+. For any e > 0 and any K <g r 3 / 2 ('w) ; there exists a radius r 0 = ro(e,K,w) £ 
(0,1), such that r 3 / 2 (w) n K is (e, ro )-Reifenberg flat. 

Proof. Step 1: Uniformity in compact sets. As in Proposition 3.5, the Corollary 
reduces to showing the following claim: 

Claim. Let aU : B+ —> be a uniformly elliptic W 1,p tensor field 

with p £ (n + l,oo]. Let w be a solution of (4) in B+. Then, for any e > 0 and 
K (g T 3 / 2 (w), there exists ro = ro{e,w,K) such that for any Xq £ K and r < ro, 
there is a matrix S = S(r,x o) £ SO(n + 1) such that 

( 14 ) IKotOSz) - w 3/2 {x) || C1(B +) < e, 
where w%/ 2 {x) is as above. 

In order to infer the claim, we argue by contradiction and assume that the 
statement were wrong. Then there were an e > 0, a sequence rfc with r*, —> 0 and 
and a sequence Xk with Xk £ K such that for all S £ SO(n + 1) 

\\ w x k ,r k (Sx) - w 3 / 2 (:r)|| C i( B +) > e for all k £ N. 

By compactness of K we may assume that Xk —> x. However, by compactness 
(which for variable centers we infer from a combination of Proposition 4.3 in 
[KRS15] and Corollary 4.1 in [KRS15]) of the sequence w Xk , rk and by Lemma 
3.3, we obtain a subsequence, which we do not relabel, and a rotation S such that 

\\w Xk ,r k {x) - w 3/2 (Sx )|| C1(B +) -l 0. 
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Step 2: Conclusion. By rescaling with the radius r o = ro(e, K,w), it is possible to 
reduce the claim of Corollary 3.7 to the setting of Corollary 3.6. □ 

3.2. Lipschitz regularity. In this section we present a first improvement of the 
regularity of the regular free boundary. 

In the sequel we will always assume that 0 G P,„. Moreover, we will use the 
notation 

w 3/2 (x) := c n Re(x„ + ix n+ i) 3 ^ 2 , 

to denote our model solution. Here c n > 0 is a normalization constant such that 
||w; 3 / 2 || i 2 ( B +) = 1. Keeping this convention in the back of our minds, we prove the 
Lipschitz regularity of the regular free boundary: 

Proposition 3.8 (Lipschitz regularity). Let a u : B+ — > be a tensor 

field, which satisfies (Al), (A2), (AS), (Af). Let w be a solution of (f) in B+. 
Assume that for some small positive constants e 3 and c» 

(i) \\w - w 3/2 \\ cl{B + ) < e 0 , 

(ii) l|Va ij '|| iP{B +) < c*. 

Then if eo and c* depending on n,p are chosen sufficiently small, there exists a 
Lipschitz function g : B"j 2 —> K such that (possibly after a rotation) 

r w n B[ /2 = {x n = g(x")\x" G By 2 }. 

Moreover, there is a large constant C = C(n,p), such that || V"g\\L°°(B'^ /2 ) < 
C maxjeo, c*}. 

We remark that in addition to the results from Proposition 3.8, we actually show 
that g - and hence the regular free boundary - is C l (B",f) regular (c.f. Remark 
3.20). 

In order to prove the Lipschitz regularity of the free boundary, we deduce positiv¬ 
ity of the directional derivatives in a cone of tangential directions (c.f. Proposition 
3.17). To achieve this, we use a comparison argument, which is standard in the 
constant coefficient case. We start by considering the equation for v := d e w, where 
e G S n fl B( is a tangential direction. From the H 2 estimate of w (c.f. [Ura87]), 
v G -/^(T?) 1 "). By differentiating the equation of w, we obtain the equation of v in 
Bi n {x n+1 > 0}: 

dia^djV = diF 1 in B+ fl {x n +i > 0}, where F l = ~(d e a l -’)djU ;. 

Compared with the constant coefficient situation, the fact that the coefficients of 
our equation are not better than W 1,p regular, leads to a divergence right hand 
side with F 1 G L p . This causes difficulties in applying the comparison argument 
directly. To resolve this difficulty we use an appropriate decomposition to split 
the solution v into an “error term” v\ (c.f. equation (17)), which deals with the 
divergence term F l , and into a “main term” v 2 (c.f. equation (18)), which captures 
the behavior of the respective solutions. 

Remark 3.9 (Reflection and extension). As it proves convenient to work on the 
whole ball instead of only on the upper half ball (for instance for applying interior 
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elliptic regularity theory), we reflect the metric a 10 to B 1 by setting 
a l] (x',x n+ 1 ) = a tJ (x',-x n+1 ), i,j = 

(15) a n+1 ’ 3 (x',x n+1 ) =-a n+1 ’ J (x',-x n+1 ), j = l,...,n, 

a n+1 ’ n+1 ( X ',x n+1 ) = a n+1 ’ n+1 (x',-x n + 1 ), 

for (x',x n +i) £ Bf. We reflect w to Bf by defining w(x',x n +i) = w(x', —x n +i). 
Here we use the same notation to denote the functions after refection. Due to the 
off-diagonal assumption (Al), the reflected coefficient matrix a 1 - 7 is in W 1,P (B i). 
Since <9„+i w = 0 in B[\ A w , we have w £ C 1,1//2 ~ s (B i \ A w ) fl H 2 [B\ \ A w ), where 
5 £ (0,1) is an arbitrarily small constant, and v = d e w satisfies 

(16) dia^djV = diF 1 in Bi \ A w , F l = —(d e a l ^)djW, 

v = 0 on A w . 

As we are mainly interested in the local properties of solutions and of the (regu¬ 
lar) free boundary and in order to avoid technical issues on the boundary dBiC\A w , 
we will often consider (16) as an equation in the whole space R n+1 : We sim¬ 
ply extend a 1 - 7 to R" +1 such that the extended coefficients remain Sobolev regular, 
a 1 - 7 g W l )(((M. n+1 ), with equivalent ellipticity constants and with || Va 1 - 7 ||lp(r«+i) at 
most twice as large as it had originally been. Moreover, we extend F l to R” +1 by 
multiplying by the characteristic function xb x ■ With a slight abuse of the notation, 
we denote the set A w D B[ by A w . In the sequel, whenever the equations are writ¬ 
ten in the whole space, we will always have this extension procedure in mind, even 
though we may not refer to it explicitly. 


In order to deal with the L p divergence form right hand side in (16), we decom¬ 
pose v = v\ + V 2 with 

dia^djVi — K dist(x, T w )~ 2 v\ = diF 1 in R n+1 \ A w , 

vi = 0 on A w , 
where AT is a large constant. Then V 2 := v — v\ solves 

dia^djV 2 = —K dist(x, T^) -2 !;! in R n+1 \ 

V 2 = 0 on A w . 


(17) 


(18) 


Intuitively, we expect that v\ is a “controlled error” (c.f. the bounds in Proposition 
3.10) and that V 2 determines the behavior of our solution v. 


That this is indeed the case and that v\ can be treated as a “controlled error” is 
a consequence of the following more general result which will be proved in Section 
3.4: 

Proposition 3.10. Let A and T be closed, non-empty sets in R” x {0} with T = dA. 
Let a 1 - 7 : R” +1 — > a bounded, measurable, uniformly elliptic tensor 

field whose eigenvalues are in [1/2,2]. For K > 0 consider the equation 

dia^djU — K dist(:r, T)~ 2 u = diF 1 + g in R n+1 \ A, 

u = 0 on A, 

with F l £ L 2 (R” +1 ), i £ {l,...,n+ 1} and dist(ar,T )<7 £ L 2 (R n+1 ). Then there 
exists a unique solution u with Vu £ L 2 (R n+1 ), dist(a;,T) _1 u £ L 2 (R ra+1 ) and 
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u(x) = 0 on A for H n a.e. x. Under the above conditions, if we additionally 
assume that for some a > 0 and p £ [n + 1, oo] 

F l dist(x,T)- a £ L p (R n+1 ), ie{l .n + 1}, 

<?dist(*,r) 1 - Ii F i - <T e L p / 2 (R” +1 ), 

then there exist c = c(n) large and Co = Co(n,p), such that for K = c(n)(a + 1), 

|it(ar)| < Co dist(x, r) 1 ” 2 p i+cr (J|F 1 dist(-, r) _<T || iP(Rn+1) 

+ ||<? dist(', r) 1 ^ _ ' <J |Lp/ 2 (R n + l ) ) . 

Moreover, u vanishes continuously as x approaches A. 

Remark 3.11. With slight modifications of the function spaces, the statement of 
Proposition 3.10 remains true for a < 0. 

Remark 3.12. If in addition £ W 1,P (B+ ) with p £ (n + 1, oo), it is possible to 
make the vanishing order of u towards A more explicit (c.f. Remark 3.27 after the 
proof of Proposition 3.10). More precisely, for p £ (n + 1, oo) we have 

K+)| < Cq dist(x, r) CT dist(x, A ) 1-2 ^ 1 (||FMist(-,r) _<T || iP(Rrl+1) 

+ llff dist(',r) 1 P <J || i p/2( R n + l)^ . 

Although we postpone the proof of Proposition 3.10 to the end of this section, 
we will in the sequel frequently apply it, for instance with A = A^, and T = r^,. 

The remainder of the proof of Proposition 3.8 is organized as follows: Working 
in the framework of general (e, 1)-Reifenberg flat sets A and T, in Section 3.2.1 
we first construct a lower barrier function (Proposition 3.16) which serves as the 
basis for the then following comparison argument. Here we deduce an (almost) 
optimal non-degeneracy condition in cones for elliptic equations with controlled 
inhomogeneities in domains with Reifenberg slit (Proposition 3.17). In Section 
3.2.2 this is then applied to prove the Lipschitz regularity of the (regular) free 
boundary for the thin obstacle problem (Proposition 3.8). 

3.2.1. Construction of a barrier function and a comparison argument. In this sec¬ 
tion we first construct a barrier function for the slit domain with Reifenberg flat slit 
(c.f. Proposition 3.16). This is crucial for the then following comparison argument 
which allows us to work with essentially critically scaling bounds. 

In the whole subsection, we work under the following assumption: 

Assumption 3.13. We consider a closed set A C K n x {0} with boundary T := 
x { 0 } A - We assume that 0 £ T and that TflBi is (e, 1 )-Reifenberg fiat. Here 
e > 0 is a fixed small constant. 

We start by considering a Whitney decomposition, {Qj}jeN, of B\ \T, i.e. {Qj} 
is a collection of dyadic cubes Qj with disjoint interiors that cover B\\T such that: 

(Wl) diam(<5j) < dist(Qj, T (~l B\) < 4diam(Qj). 

(W2) If Q\ and Q 2 touch, then (l/4)diam(Q 2 ) < diam(Qi) < 4diam(Q 2 )- 
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(W3) There are at most (12)" +1 cubes in {Qj} which touch a given cube Q £ 
{Qj} (cf. Chapter VI [Ste70] for the existence of such cubes). 

Moreover, as A, T C R" x {0}, it is further possible to choose the Whitney cubes 
symmetric with respect to the plane {x n +i = 0}. 

For each Whitney cube Qj with center x 3 and diameter rj , let x 3 £ T be a (not 
necessarily unique) point such that dist(xj, T) = \xj—Xj\. Using this point, we asso¬ 
ciate a rotation Sj to each Whitney cube Qj by choosing one of the (not necessarily 
unique) rotations 6'(64r j , Xj ) which are defined by the property that S(64rj,Xj)e n 
determines the normal of the plane L(64rj,Xj) in the Reifenberg condition (7) with 
S = e for some t<Cl. We define Vj := Sje n , which we refer to as the approximate 
normal at Xj in B^Axj) associated with Qj. Note that Vj £ S n D {x n + 1 = 0}, as 
we interpret T as a subset of 1" x {0}. 

We observe that the approximate normals only vary slowly on neighboring Whit¬ 
ney cubes: 

Lemma 3.14. Let A,T be as in Assumption 3.13. Let {Qj} be a Whitney decom¬ 
position of B\ \ r and let Qj, Qk be neighboring Whitney cubes. Let Vj and Vk be 
the associated approximate normals defined above. Then 

|Vj - v k \ < Ce. 

Proof. Let Xj and Xk be the centers of the Whitney cubes Qj and Qk respectively, 
and let Xj , Xk £ T,„ be (not necessarily unique) points which realize the distance 
of Xj,Xk to r w . Let rj be the diameter of Qj. By using (Wl) and (W2), it is 
immediate to check that 

0 < |Xj — Xk\ < 32 rj. 

Let Sj = S(64rj,Xj) and let v 3 = Sje n be as above. Let Lj = {x £ R" x {0}| x ■ 
Vj = 0} be an (n — l)-dimensional hyperplane in R" x {0}. Then condition (7) in 
Definition 3.4 gives that 

dist// f I',,, D - {x 3 ), Xj -I- Lj n B( 34 r^ ( x :i )) V f'ery. 

Similarly at Xk we have 

dist// (T^ n (xk )t Xk Lk H B^ rk (x/c)) V Cerk • 

Comparing the above two inequalities and using the triangle inequality for distjj, 
we deduce 

Wj - Vk\ < Ce. 

□ 

With this auxiliary result at hand, we construct a first barrier function for the 
divergence form operator L = dia^dj with a 1 - 7 satisfying the assumptions (A2), 
(A3), (A4). The idea is first to locally construct subsolutions in each Whitney cube 
Qj and then to patch them together by a partition of unity. The patched function 
remains a subsolution as the approximate normals vary slowly due to Lemma 3.14. 
In order to implement this idea, in a first step, we view the operator L as a pertur¬ 
bation of the non-divergence form operator Lq := a lJ dij. Then in the second step 
we use the splitting technique from Proposition 3.10 to derive a barrier for the full 
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operator L. 

In order to simplify notation, we introduce the following parameters as abbrevi¬ 
ations: 

Notation 3.15. We abbreviate by setting 


• £ 0 := (2 i y /n) h 


Using this notation we proceed with the barrier construction. 

Proposition 3.16 (Barrier function). Let A,T,e be as in Assumption 3.13 and 
let c*, be as in Notation 3.15. Then for any s £ (0,1/2), if e = e(n,s) and 
c* = c*(n,p, s) are chosen sufficiently small, there exists a function hj : B\ — > M 
such that 

(i) Lhj(x) > c s dist(x, r) _ i + i in B\ \ A for some c s > 0, 

(ii) hj(x) > c„ dist(x, r)5+f onB\C\{x\ dist(x, A) > £ 0 dist(x, T)}, 

3 n +1 

(iii) hj(x) > — c n dist(x, T) 2 * in B i and hj(x) =0 on A. 

Proof. Step 1: Barrier for Lq. Let rjj be a partition of unity associated to the 
Whitney decomposition of £?i\r from above, i.e. the functions rjj satisfy 0 < rjj < 1, 
rjj = 1 on \Qj and r\j = 0 for all x £ Qk if Qk £ Af(Qj) '■= {Qi\ Qi touches Qj}. 
We consider 



where v7 is constructed as follows: We first define an affine transformation, Tk, 
associated with each Whitney cube Qk- More precisely, we set T k (x) := RkBf 1 (x — 
Xk), where as before Xk denotes a point on r D B\ which realizes the distance of 
the center x k of Qk to F D B i, BkB^, := {a i: >{xk)), and Rk £ SO(n + 1) is defined 
such that RkB^Uk = Ci,k e n and RkBj.e n +i = C 2 ,k&n+ 1 - Here v k is the associated 
approximate normal which was defined before and described in Lemma 3.14. The 
constants Ci.k-, c- 2 ,k with Ci.k-, C 2 ,fc £ (-^, y/2) only depend on the ellipticity constants 
of a l i. Next, for s £ (0,1) we consider 

V (x) .— Hh/2 (^n) •En+1 ) 5 1Cl/2 (*^n 7 ) — He(x rl T i |x n _|_i |) ^ . 

A direct computation leads to 


Av (x) = s(l + s)(x 2 n +x 2 n+1 ) 1/2 (wi/ 2 ) 1+8 
> s(l + s)(x 2 n +x 2 n+1 )-i + i. 


Then we define 


v k (x) ■= v (T k x). 


We observe that by definition 

(20) Tfc(x) • e n = c“J.(t - x k ) ■ v k , T k (x) ■ e n+1 = cfl(x - x k ) ■ e n+1 . 
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Thus, v k (x) = (Re(c 1 j,(a; - Xk) ■ v k + i\c 2 j,(x - Xk) ■ e n +i|) * ) 1+s - By using the 
computation for Av~(x), one can check that v k satisfies 

a lJ (x k )dijV k (x) > s(l + s)((T k (x) ■ e n ) 2 + (T k (x) ■ e n+ i) 2 )“3+f 
in R n+1 \ {x n+ \ = 0, {x - x k ) ■ v k < 0}, 
v k (x) = 0 on {x n+1 = 0, (x - x k ) ■ v k < 0}. 

Due to the (e, 1)-Reifenberg condition on T, the definition (20) and due to the 
property (W2), there exists an absolute constant c £ (0,1) such that 

a 13 {x k )dijV k (x) > cs(l + s)r k 5+!? , r k = diam(Q fc ) 

for each x £ supp(? 7 fc). By (Wl) and (W2) this can further be rewritten as 

(21) a lj (xk)dijv k (x) > cs(l + s) dist(x, T) _ 3+f 5 x e supp ( r q k )- 

We compute Lo(Efc r !k v k )i which is 


L o(%2 VkV k ) = ^2a Z3 {x k )d ij (r]kV k ) + ^](a u (x) - a lJ (x fe ))%(r? fe x fc ) 

k k k 

= J2( alJ ( x k)dijV k )r] k + ^^(a lJ [x k )dijr]k)v k 
k k 

+ 2 y a lj ( Xk)djT]kdjV k + y^(a v (x) - a lj (x k ))dij (rj k v k ). 
k k 

Given any x £ B\ \ T, there is a unique Qe such that x £ Qe- By definition of 
our partition of unity, rjj ^ 0 in Qe iff Qj £ AT(Qe). Using (21) and (W2) we thus 
deduce (with a possibly different c) 

y^(q ZJ (xk)djjV^ ( x))rjk (x) > cs( 1 + s)r^ 1+ ^. 
k 


This will be the main contribution in L^h . In order to see this, we now estimate 
the error terms: For any x £ Qi : 

(xk)dijTik(x))v k (x) 
k 

= 'Yy? 3 {x k ) - a %3 {x i )\d ij r] k (x)v k {x) + ^a l3 {xi)d ij r}k(x)v k {x) 
k k 

= ^[a u (x fc ) - a 13 (x e )]dijrik(x)v k (x) + 'SC a 13 ( x e )d ij r] k (x) (x) - vj(x)\ 


<C(n) (, 


-y———H - 

c * r i 2 2 + er, 9 


3 

2 T 2 


7 = 1- 


n +1 
P 


where in the second equality we used a 13 {xi)dij(J2 k Vk(x)) = Efc ^( x e)dijrjk{x) = 
0. In the third inequality we applied Lemma 3.14 and noted that for x £ supp(%) 
and Q k £ Af(Qe), 


(22) \v k (x) - v e (x)| = \v (Tfe(x)) - v (7)p(x))| < C max{ c»rj, e}r e 


Indeed, this follows, since |7i-(x) — T;(x)| < max{c*rj, e}r^. 

Similarly, by using Efe a *' ? ( x t)^ 7 7fc( a: ) = 0 for &U j 6 (l,...,n+ 1} as well as 
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Lemma 3.14, we infer that 

( x k)diVkdjVk 

k 

= - a lJ {xe)\diVk{x)djV^ (x) + ^ a %J {xe)d i rik(x)[d j v^ (x) - d-jvj (x)] 

k k 

< C(n ) (c*r; -i+f +er7 i+f ) . 

Analogous to the estimate from before we used that 

s — 1 

I dj v k( x ) — djvj{x )| = | djV~(Tkx) — djV~(T^x)\ < max{c»rj, e}r e 2 
in the last inequality. Finally, 

(x) - a lJ (x k ))dij ( r) k v £ (x)) < C(n)c»rJ 2+ 5 . 
k 

Combining all the previous estimates, if e = e(n,s) and c* = c»(n, s) are small 
enough, we have for some c s > 0 

(23) L 0 h~(x) > c s dist(x,r)“5+2. 

Next we show that h~ = 0 on A. By construction we directly note that h~ = 0 
on A \ r. Due to the continuity of hr and due to the fact that F = 9A, this then 
also holds on T. 

Step 2: Perturbation. We now modify h ~, so as to become a barrier for the full 
operator L. For this we note that G = XB 1 (dia i:i )djh~ satisfies the bound 

||G(x)dist(x,r)5|| Z/ p (R n + i) < C(n)c», 

since 

|<9,7i“(x)| < C(n) dist(x, r) _ 5 f or a n x £ B\ \ A. 

Let q(x) be the solution of 

Lq{x) — K dist(x, r)~ 2 < 7 (x) = — G(x ) in R" +1 \ A, 
q(x) = 0 on A, 

for some large constant K = K(n). Then, by Proposition 3.10 there exists a 
constant Co(n,p) such that 

(24) k(x)| < Cbc* dist(x, T) 
and q vanishes continuously up to A. Let 

h~{x) := h~(x) +q(x). 

Thus, for c* = c*(n,p, s) small, we have 

L(h~) = L 0 h~ + {dia l ^)djh~ + ATdist(x, P)~ 2 q — ( dia z;i )djh~ 

> c s dist(x, F)~^ + i — CqKc* dist(x, r) _3_ — 

> ic s dist(x,r)^^ + 5. 

Moreover, we note that for each k with supp(? 7 fc) fl {dist(x, A) > tb dist(x, T)} ^ 0 
and for each x £ supp(ryfc) fl (dist(x,A) > £q dist(x, T)}, we have that (x) > 
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c n dist(x,F) 2 + 2 . Using the definition of h and the estimate (24) for q(x) (as well 
as a sufficiently small choice of c*), we therefore arrive at 

(25) h~(x) > c n dist(ai,r)5 + 2 for x £ B\ n {dist(x, A) > £q dist(x, T)}. 

□ 


Using the previously constructed barrier function, we proceed by deducing a 
central non-degeneracy condition (in cones) for solutions of elliptic equations with 
controlled inhomogeneities in domains with Reifenberg slits. This corresponds to a 
comparison principle for solutions to divergence form, elliptic equations with rough 
metrics and divergence form right hand sides in domains with Reifenberg slits. 

For convenience of notation, we change the geometry slightly, i.e. in the following 
we use cylinders instead of balls: For r, £ > 0 and {0}, B' r (xo) x {—£, £) := 

{{x',x n+1 ) £ R n+1 | \x' - x' 0 \ < r, |z„ + i| < £}. 


Proposition 3.17. Let A, T, e be as in Assumption 3.13. Assume that c„, £q are 
as in Notation 3.15. Suppose that u £ H 1 (B[ x (—1,1)) n C(B[ x (—1,1)) solves 

dia^djU = diF 1 + gi + g 2 in B[ x (—1,1) \ A, u = 0 on A, 


and that the following conditions are satisfied: 

(i) The metric a 13 satisfies the assumptions (Al), (A2), (AS), (A)) and for 
some constants 5 q £ (0,1), a > \ and 6 > 0, 


n+1 


Y ||dist(-,r)- CT F l | 


L*(B(X(-1,1)) 


dist(-,r ) 1 ? a gi 
< 6 . 


LP/ 2 {B[x(- 1 , 1 )) 


dist (•, T) 2 5 °g 2 

L°° {B[ X (— 1 , 1 )) 

(ii) u > 1 on B[ x [—1,1] D {|o;„+i| > £q}. 

(iii) u > —2~ 8 in B( x (—£ 0 ,£ 0 ). 

Then, if e, S, c» are sufficiently small depending on n,p, So, a, there exists a constant 
c n > 0 such that 


(26) 


u(x) > c n dist(x, r) 2 + 2 , e = min{(5o, ^ + cr} 

2 p 


for all x € {B(^ 2 x (—1/2,1/2)) fl {a:| dist(x, A) > £ 0 dist(x, T)}, and 


(27) u(x) > — c„ dist(x, T) 1 p +<T 

for all x £ B( /2 x (-1/2,1/2). 

Remark 3.18. By a scaling argument, Proposition 3.17 remains true, if condition 
(iii) is replaced by u > cq on B( x [—1,1] fl {|x n +i| > £o} for some cq > 0 (and if 
the remaining constants in Proposition 3.17 are rescaled appropriately). 


Remark 3.19. The parameter a > — 1 quantifies the necessary decay con¬ 

ditions on the inhomogeneities F l , i £ {1,..., n + 1}, g±, g 2 in dependence on 
the integrability of these quantities (measured in terms of p). In the application 
to the thin obstacle problem, the corresponding inhomogeneities satisfy these decay 
assumptions, if either 

• p > n + 1 and no obstacle is present, 

• or p > 2(n+ 1) in which case an arbitrary W 2,p obstacle is allowed. 
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More precisely, in the setting of the thin obstacle problem, we apply Proposition 3.17 
to the tangential derivatives d e w, which satisfy the divergence form equation 

dia^djU = diF 1 , where F z = (d e a^)djW + d e (a^djcf), 

with metric a lJ £ W 1,p and obstacle (f £ W 2,p . We note that on the one hand by 
Lemma f.3 in [KRS15], the first term in F l satisfies 

dist(ir, T w )~ 1 ! 2 ln(dist(:r, T w ))~ 2 {d e a^)djW £ L p . 

On the other hand, the second term d e (a lJ dj(f>) £ L p , which originates from a 
potentially non-vanishing obstacle (c.f. Section 5.2), does not carry any additional 
decay. However, 

• if p £ (n+ 1,2 (n+ 1)], then a > 1 > 0, i.e. decay on F 1 is required. 

The assumption (i) from Proposition 3.17 on F l is thus in general only 
satisfied if we set (f> = 0. 

• If p £ (2 (n + l),oo], we have ^ < 0. Hence, we do not need any 

additional decay on F 1 . In this case Proposition 3.17 applies for any <fi £ 
W 2 ’ p . 

Proof. In order to prove the result, we consider an appropriate comparison function. 
For xq £ 11^2 x (—£o,£o), we set s = e (where e is the constant from (26)) and 
consider 

h(x) := u(x ) + \x' - (x 0 )'\ 2 - 2~ 8 hf(x) - 2 nx 2 n+1 . 

A direct computation shows that h satisfies 

dia^djh = gi+g 2 , 

where 

gi = diF 1 + gi+ 2 (xj - {x 0 )j)(dia lJ ) - 4n(x n+ i)(5 J a ra+1 ’- 7 ), 
g 2 = g 2 - 2 - s dia ij djhf + 2a" - 4na n+1 ’ n+1 . 

We split h = hi + h 2 , where 

dia^djhi — K dist(:r, T)~ 2 hi = gi in R” +1 \ A, 

hi = 0 on A. 

By Proposition 3.10, there exists Co = Co{n,p) such that 
(28) \hi(x)\ < Co(c* + 5) dist(:r, r) 1- ~ hCT . 

Now we consider h 2 . Recalling the definition of £q, we have that 

h 2 > 1 — C 0 (c» + 5) — 2~ 8 — 2 nil > ^ on B' 3 x {±4}, 

2 4 

h 2 > —2~ 8 + - C 0 (c» + 5) — 2~ 8 - 2u£q > 2~ 8 on SB's x (-4,4)- 

lo 4 
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We further recall that u = hi = h~ = 0 on A. Thus, h 2 > 0 on A. 

If 8 and c* are chosen sufficiently small, h 2 satisfies 

dia^djhi = g 2 — K dist(x, T)~ 2 hi 

< {8 + c*)CqK dist(x, r)^ 1_ ~ +<T — 2~ 8 c s dist(x, r) _ 2 + 2 

+ 2a u - 4na n+1,n+1 + 8 dist(x, T)~ i +<5 ° 

<0 in (Bi x (—4,4))\A, 

4 

where we used the smallness of 8 and the definition of s in terms of e. Hence, by the 
comparison principle, /12 > 0 in B'3 x (— £q,£o)- Combining this with the estimate 

(28) gives 

h{x) > —Co(c* + (5) dist(x, r) 1_ ~ +;T in B'3 x (—^o,^o)- 

4 

In particular, 

h(x 0) > — Co(c» + 8) dist(xo, T) 1 ” - ? h<T . 

Rewriting this in terms of u 2l yields 

u(x 0 ) > -C 0 (c* +8) dist(x 0 ,T) 1 %- +a + 2~ s hj(x 0 ). 

Using the growth properties of the barrier function h~ (c.f. Proposition 3.16) 
and by choosing c*,(5 sufficiently small in dependence of n,p and <5o, implies the 
estimates (26) and (27) in B' x / 2 x (— £q,£o). Finally, the full non-degeneracy in 
B' x j 2 x (—1/2,1/2) fl {x\ dist(x, A) > dist(x, T)} follows from our assumption (ii) 
in combination with the fact that 0 £ T. □ 

3.2.2. Proof of Proposition 3.8. With the aid of Proposition 3.17 we can now prove 
the Lipschitz (and C 1 ) regularity of the regular set of the free boundary. We 
only sketch the proof as there are no real modifications with respect to the one in 
[PSU12]: 

Proof of Proposition 3.8. By choosing eo and c* sufficiently small and by invoking 
Proposition 3.5, we may assume that P^ fl B[ , 2 is (e, 1/2) Reifenberg flat. Using 
this and the C X (B X ) closeness of w to w 3 / 2 , we apply Proposition 3.17. In order 
to satisfy its assumptions, we transfer positivity from d e w 3 / 2 to d e w, where e £ 
S n ~ 1 x {0} is a tangential direction. We consider £0 = (2 4 v / n) _1 as in Proposition 
3.17. Then for any 77 > 0 and any tangential vector e £ {x £ S n fl B[ | x n > r)\x"\\, 
the function w 3 / 2 satisfies: 

d e w 3 / 2 (x) > 0 if x £ Bi, 

d e w 3/2 (x',x n+1 ) > c n r](l + rf‘)~ l ^ 2 £ 1 J 2 > 0 if x £ B x n{|x„+i| > £ 0 }. 
Moreover, d e w solves 

di(a 1 ^djd e w) = diF 1 in Bi \ A w , F l = ~{d e a l ^)djW. 

By Lemma 4.3 in [KRS15], || dist(-, Tu,) - ^ ln(dist(x, r ? „)) _2 F i || i p( Bl ) < C(n,p)c*. 
Hence, by Proposition 3.17 (applied with <7 = \ — /x, where p is an arbitrarily small 
but fixed constant), if eo and c* are chosen sufficiently small depending on n,p, p, 
then there exists a positive constant c RiP such that 

(29) d e w(x) > c n , P dist(x,ru,) 5+ 2 > 0 
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for all x £ B 1 / 2 C\{x\ dist(x, A™) > £q dist(x, P„,)} and e € {x € S n (lB' 1 \x n > rj\x"\}. 
Here 17 > Cmax{eo,c*} for a large enough positive constant C = C(n,p) and 
e = a + \ — Moreover, Proposition 3.17 also implies the global lower bound 

d e w(x ) > — dist(x, T ^,) 1 p h<T in B 1 / 2 ■ As d e w = 0 on A w and as fl w is covered 
by the set (x| dist(x, A w ) > £0 dist(x,F UJ )} (on which the bound (29) holds), we 
deduce that d e w > 0 in with strict positivity in By, 2 D {x| dist(x, A„,) > 

£0 dist(x, Tu,)}. As a consequence (c.f. [PSU12]), 

{w(x r , 0) > 0}nBj /2 = {x n > g(x") | x" £ B" /2 }, 

where g is a Lipschitz continuous function with \\A7"g\\L°°{B” ) < C max{ eo, c* }. 

1/2 

□ 

Remark 3.20. If we return to an arbitrary solution w of (4), we can even con¬ 
clude the C 1 regularity of the free boundary from the previously given argument. 
Indeed, after rescaling, it is always possible to satisfy the closeness assumption (i) 
in Proposition 3.8 with an arbitrarily small parameter e$. Hence, we can infer that 
for each 77 > 0 there exists p v > 0 such that ||V ,, 5 || ioo ^ B „ ^ < p. This yields that g 

is differentiable at 0 , or in other words, the tangent plane of T^i 2 {w) exists at the 
origin. The same holds true at other regular free boundary points. Moreover, from 
(29) it is not hard to see that the tangent plane varies continuously in a neighbor¬ 
hood of the origin. Hence, we indeed obtain that r 3 / 2 (ir) C\B p is not only Lipschitz 
but C 1 regular. 

Remark 3.21. At this stage and with these (comparison) techniques we do not 
know how to further estimate the modulus of continuity. Heuristically, the im¬ 
provement comes from the fact that the convergence to the homogeneous solution is 
much faster than the compactness arguments suggest, once the solution is close to 
the homogeneous solution in the sense of Proposition 3.5. Thus, in the next section 
we use a boundary Harnack inequality to overcome this difficulty. 

3.3. C 1,a regularity. In this section we improve the regularity of the (regular) 
free boundary to obtain its C 1,ot regularity for some a £ (0,1]: 

Proposition 3.22. Let be a uniformly elliptic W 1,P (B((), p £ (n + l,oo], tensor 
field with a l:, (0) = S lJ . Assume that the conditions (i) and (ii) of Proposition 3.8 
are satisfied with the same constants eo and c*. Then if eo,c* are sufficiently small 
depending on n,p, there exists a = a(n,p) £ (0,1 — such that g £ C 1 , “(B"y 4 ). 
Moreover, ||V"g|| c o,°(B" 4 ) < C(n,p) max{e 0 , c*}. 

We remark that Theorem 1 is an immediate consequence of Proposition 3.22: 

Proof that Proposition 3.22 implies Theorem 1. As a result of the analysis in [ vRS15], 
we may assume that for an appropriate sequence of radii Vj > 0 , rj —> 0 , the 

(L 2 )-rescalings, w rk {x) := n+ i"( rfc —-, converge to a global solution wo in 

r k 2 IMI LHB+ k ) 

c'L“( r + +1 ) for some a £ ( 0 , 1 / 2 ), where after a possible rotation of coordinates 
w 0 (x) = w 3 / 2 (x) = c n Re(x n +ix n+1 )T 
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For each k , w rk is a solution of (4) in B r -1 with metric a]? , where a l f k (x) = a 1 - 7 (r^x), 
with rescaled free boundary r„,. A direct calculation shows that 

l|Vo«|U, (fll) =ri" 2 ^||Va«|| iP ( BrJ . 

Thus, for rfe 0 small enough, assumptions (i), (ii) in Proposition 3.8 are satisfied for 
w rk . Thus, Proposition 3.8 implies that P uv D B l / 2 is a Lipschitz graph with 
sufficiently small Lipschitz constant. Hence, Proposition 3.22 is applicable to w rkg . 
Rescaling back to w yields Theorem 1. □ 

The improvement from Lipschitz to C 1,a regularity of the regular free boundary 
is achieved with the aid of a boundary Harnack inequality (c.f. Theorem 7 in 
[ACS08]) applied to tangential derivatives v = d e w in the slit domain Bi\A w . Here 
we follow the ideas of [CSS08] by making use of the non-degeneracy condition of the 
solution. In this context, we have to be slightly more careful, as our inhomogeneity 
is in divergence form. Consequently, we rely on the decomposition, v = iq + V 2 , 
which we introduced above (c.f. Proposition 3.10). Another difficulty, which we 
have to overcome, is the fact that we may not assume that v is positive in the whole 
domain B\ \ A w : On the one hand it satisfies (by Proposition 3.17) 

v(x) > c 0 dist(a:,r iu )3 + * in {x € £?i| dist(x, A w ) > £ 0 dist(x, r„,)}, 

for £o = (2 4 v / n) -1 and e = min{<5o, cr + i — while on the other hand in the 

complement we only have 

v(x) > —ci dist(x,r w ) 1 ~~p- +r7 in B\. 

Here a > — i is the same constant as in Proposition 3.17 (i). 

Despite these additional difficulties, it is possible to adapt the general strategy 
of [CSS08]: Working in the framework of general elliptic equations with controlled 
inhomogeneities in domains with Lipschitz slits, in Section 3.3.1 we first prove a 
Carleson estimate (Lemma 3.24) which is adapted to our setting. Then we continue 
with an appropriate boundary Harnack inequality (Lemma 3.25). In Section 3.3.2 
we finally apply these results in the setting of the thin obstacle problem and hence 
infer the desired regularity result of Proposition 3.22. 

3.3.1. Carleson and boundary Harnack estimates. In the sequel, we work in the 
slightly more general set-up of elliptic equations with controlled inhomogeneities in 
domains with Lipschitz slits. Hence, in this section we make the following assump¬ 
tions: 

Assumption 3.23. Let g : R” -1 —> R be a Lipschitz function with Lipschitz 
constant e > 0 and g( 0) = 0. Then we consider 

A := {x n < g{x")} Cl"x {0}, 
r := {x n = g(x")} Cl"x {0}. 

We remark that the Lipschitz regularity of T in particular implies its (e, 1)- 
Reifenberg flatness. Moreover, in our set-up, the smallness assumption on the 
Lipschitz constant of g does not pose restrictions and can always be achieved by 
choosing the constants eo, c* sufficiently small in Proposition 3.8. 
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As in the previous Section 3, we use the abbreviations from Notation 3.15, i.e. 
we again set Iq := (2 4 - v /n)~ 1 and c» := ||Va lJ ILp(_b 1 ) f° r notational convenience. 

In this framework we derive the following Carleson estimate: 

Lemma 3.24 (Carleson estimate). Let A, T, e satisfy Assumption 3.23. Suppose 
that u £ C(B i) D H 1 (Bi \ A) solves 

dia^djU = diF 1 + gi + g 2 in Bi \ A, u = 0 on A, 

where we assume that condition (i) from Proposition 3.17 and the following non¬ 
degeneracy assumption holds: For a , e as in Proposition 3.17 we have 

u(x) > dist(x, r)5 + S in {x £ BA dist(a?, A) > £q dist(a;, T)}, 

(30) „ +1 

u{x) > — dist(x,r ) 1 p +<T in Bi. 

Then if 5 = S(n,p,5 0 ) and e = e(n,p,5o) are small enough, there exists C = 
C(n,p,e) such that 

sup u < Cu(Q + re n ), 

BAQ) 

for all Q £ r (~l B ^ 2 and r £ (0,1/4). 

Proof. Let Q £ T D B^j 2 and r £ (0,1/4). Let u* be the solution of 

dia^djU* = 0 in B 2r (Q) \ A, 

u* = max{«, 0} on d(B 2 r (Q) \ A). 

Then on B r (Q) \ A Theorem 7 of [ACS08] implies that there exists C = C(n,e) 
such that 


u*(x) < Cu*{Q + re n ). 

We transfer this estimate to u by exploiting a comparison argument which makes 
use of the non-degeneracy of u. We consider 

h(x) := u*{x) — u(x) — h~{x ) + 8 r 5 + 2 , 

where h~ is the barrier function from Proposition 3.16 and s = e. For h we obtain 
dia lj djh = —diF -1 - gi - g 2 - dia lj djh~ in B 2r {Q ) \ A, 
h > 4r (1+s)/2 on d(B 2r (Q) \ A), 

where in the second inequality we have used \u* — u\ = max{— u, 0 } < r p +a on 
d(B 2r (Q) \ A). In order to estimate u, we use the splitting technique again. With 
a slight abuse of notation we write A = A fi B 2r (Q) for simplicity. We decompose 
h = hi + h 2 , where hi solves 

dia^djhi — K dist(a;, V)~ 2 h\ = — diF 1 — g\ in B\ \ A, 

h\ = 0 on A. 


Then by Proposition 3.10, 

|/ii(x)| < C(n,p)5dist(x, T) 1 p l_<J . 
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For 5 = S(n,p) sufficiently small, h 2 > 0 on d(B 2r (Q) \ A). In B 2r (Q) \ A, for 
5 = S(n,p,So) sufficiently small, 

dia'tdjhn = —K dist(x, V)~ 2 hi — g 2 — djO^ djhf 

< KC(n,p)S dist(x,T)~ 1 ^^ 1-17 + S dist(x, T) _ * +<5 ° 

— c(e) dist(ir, F) _ 5 +e 

< 0 . 

Thus, by the comparison principle we have h 2 > 0 in B 2r (Q). Hence, we obtain 
h{x) > —r 1__ p _+<T for any x £ B 2r (Q) \ A. 

Combining this with the choice of s yields 

u* — u> —Cr^ + i in B 2r (Q) \ A, 

for some absolute constant C. Similarly, we deduce u — u* > —Cr^ + % in B 2r (Q)\A. 
As a result, for x £ B r (Q), 

u{x) < u*(x) + Cr^+t 

< Cu*{Q + re n )+Cr?+i 

< Cu(Q + re n ) + Cr? + % + CV* + * 

< 3 Cu(Q + re n ), 

where in the last line we used the non-degeneracy of u. □ 

With this preparation, we come to the boundary Harnack inequality. Let A and 
T be as in Assumption 3.23. We now in addition choose 0 = 9(n,e) £ (0, 27r) such 
that for any Xq £ T 

xo +C' s (e„) C {x £ 1" x {0}| dist(x, A) > £o dist(ir, T)}. 

Here C' g (e n ) is a cone in l n x {0} with opening angle 9 and £q is as in Notation 
3.15. 


Lemma 3.25 (Boundary Harnack). Let A, T, e satisfy Assumption 3.23. Suppose 
that ui,u 2 £ C{B\) fl H l (B\ \ A) solve 

dia^djUi = diF 1 + gi + g 2 in B i \ A, u\ = 0 on A, 

dia i: *djU 2 = diF 1 + 51+32 in Bi\A, u 2 = 0 on A, 


where 

• the coefficients a lJ and inhomogeneities satisfy the condition (i) from Propo¬ 
sition 3.17 with the same constants a and S, 

• both functions u\, u 2 satisfy the non-degeneracy condition (30) of Lemma 3.2f 
with the same constant e. 

Then, if e, 5 and c* are sufficiently small depending on n,p and do, there exists a 
constant s 0 = So(n,p, e, a) > 0, such that 


so 


n 2 { 2 ^n) 


u 2 (x) _ l u 2 (\e n ) 

u\(x) ~ 0 ui{\e n ) 


in 1J (x 0 +C' e (e n )) (~1 H 1/2 . 

a: o^-6i/2nr 
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Moreover, U 2 /U 1 extends to a function which is Holder continuous on r D 
More precisely, there exist constants a G (0,1] and C > 0 depending on n,p and e, 
<7 such that for all xq G T n H(y 4 


u 2 , s u 2 , x 

— W- i x o) 

U\ U\ 


< c M\e n ) 

Ul( ~2^" n ) 


x - x 0 \ , 


x € B 1 / 4 (x 0 ) fl (so + C' e (e n )). 


Proof. We only prove the Holder continuity at so = 0. The estimate at the remain¬ 
ing points follows analogously. 


Step 1: Boundedness of the quotient. Without loss of generality we may assume 
that ui(^) = u 2 ( ~^ 2 ~) = 1. Then the Carleson estimate (Lemma 3.24) and the 
non-degeneracy assumptions on ui,u 2 imply that for all 5 > 0 

ui(x) < C in B 1 , 

u 2 (x) > (.q + 2 for all s with |s n +i| > 4 >- 

Therefore, Proposition 3.17 implies that there exists sq € (0,1) (depending only on 
the admissible quantities) such that 

u 2 — SqUi > 0 in {s G Bi | dist(s, A) > £q dist(s, T)}. 

Hence, in particular, 

— 7—7 > so for all x G {s G Si | dist(s, A) > £ 0 dist(s,T)}. 
u i(s) 2 

Exchanging the role of u\ and u 2 , we have 


Ui — SqU 2 > 0 in {s G Hi | dist(s, A) > £q dist(s, T)}. 


Step 2: Iteration argument. We only prove the Holder regularity at so = 0. 

We show that there exist a k , b k depending on s 0 and n,p, e such that 

(a) s 0 < a k < 1 < b k < s ^ 1 and b k — a k < Cp\ for some G (2~%, 1), 

(/3) there exists p 2 G [ 2 _ 2 ; p, 1 ] with b k — a k > Cp 2 for some constant C > 4, 
( 7 ) on B 2 -k nCg(e„) we have a k u\ < u 2 < b k u\. 

We argue by induction and define 

_ , x u 2 (2~ k x) - a k ui(2~ k x) _ b k u\{2~ k x) — u 2 (2~ k x) 

w i(s) :=- - -, w 2 (x) := 


b k a k 


b k - a k 


As ihi(s) + u> 2 (x) = u i(2 k x ), we may without loss of generality assume that 
wi(^f) > \u\{ 2 2 e " ). Dividing by ui( 2 2 e " ) and setting 

U\(2~ k x) 


x w i(s) 

w i (s) := , 2 „ fcen 


U\{- 


y 


i(s) := 


Ul( 


2 ~ k e 


0 


we obtain that 


2w 1 (^)>u(^) = l. 


Seeking to apply step 1, we have to show that 2uq and u satisfy its assumptions, 
i.e. the bounds on the associated inhomogeneity and the non-degeneracy. For this 
we observe: 
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• The inhomogeneity associated with 2wi satisfies the bounds from Lemma 
3.24. Indeed, w\ solves 


dj(ilj_ k djVJi = diH 1 + G\ +G 2 in Bi \ A 2 -* *, 


where 


H\x) = 


Gr(x) = 


G 2 (x) = 


>-k 


(■bk ~ a k )ui(2-^) L 


F\2~ k x) - a k F\2~ k x) 


-k„ 


>-2fc 


(b k -a k )u i(^) 


■)— 2 /e 


(ftfc - 


[51(2 fc a;)-a fc gi(2 fc x)] 

[g 2 (2~ k x) - a k g 2 (2~ k x)\ 


0-2-k(x) = a lJ {2 k x ), 

A 2 -* = 7 ^; = { x G fi il 2 ” fe;r € A i’ 
r 2 -* = ^- = {ze£M 2 _fc x e r}. 

Since dist(x, T 2 -fc) = 2 k dist(2 _fc x, T), we estimate 


(31) 

[|dist(-,r 2 - fc )- CT u i | 

2 -fe(l-2±i+a) 


LP(B 1 ) 


dist(-,r 2 - fc ) 1 - —- CT Gi 


LP/ 2 (B 1 ) 


dist(-, r 2 -fc) 2 s °g 2 


< 


{b k -a k )u i(^) 


dist(-,T) < 7 F l || i p (s )+a fe ||dist(-,r) F l \\ LP{B ) 


+ 11 dist(-,T ) 1 p a gi\\Lp/ 2 (B 2 _ k ) + °fcll dist(-,T ) 1 p a gi\\Lp/ 2 (B 2 _ k ) S ) 

2~k(^+So) 


(b k -a k )u i(^) 


(||dist(-,r)f So h\\L°°(B 2 _ k ) +afc||dist(-,r)f 5 ° g 2 \\ L °°(B 2 _ k )) 


< 


(b k -a k )u i(^) 


-25- 2 “ fc( 5 + ? ) < , 5 . 


Here we used the non-degeneracy assumption (30) in the form ui ( 2 2 en ) > 
2 -( fc + 1 )(l+f) (for which we recall that e = min{cr + i — 2 +i,(io}) as well 

as the condition (/?), i.e. £12 > 2 _ i. The inhomogeneity associated with u 
can be estimated similarly. 

• 2wi and u satisfy the non-degeneracy assumption. As 2w\{e n /2) > 1 and 
u(e n /2) = 1, Harnack’s inequality in the domain B\ D {x| dist(x, A 2 -ic) > 
£0 dist(x, T 2 -fe)} combined with Proposition 3.17 (after a rescaling which 
only depends on n) implies that there exists Co = co(n,p) such that the 
non-degeneracy condition (30) is satisfied for 2w\ and u: 


2w\(x),u(x) > Codist(x,T 2 -fe ) 2 + 2 

in {x € Bi\ dist(x, A 2 -k) > £0 dist(x, T 2 —#=)}, 

2wi(x),u(x) > —Co dist(x, T 2 -fc ) 1 p~ +a 

in {i£ Bi\ dist(x, A 2 -&) < £ 0 dist(x, T 2 -fc)}. 


L°°(B 1 ) 
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We remark that the constant Co can be chosen in the same way in each 
iteration step. 

Therefore, step 1 implies (with a possibly different so but with the same parameter 
dependence) 

> so in {x £ £?i/ 2 | dist(x,A 2 -fc) > dist(x, T 2 -ic)}. 

Rescaling back and using the scale invariance of C' e (e n ), we infer 


Therefore, 


u 2 - a k u i 
(bk - a k )u\ 


> y in B 2 -k -1 n C' e (e n ). 


(a k + y (bk - a k ))u! < u 2 < b k ui in B 2 - k -i n C e {e n ). 

Thus, we define afc+i := a k + 4y( b k — a k ) and b k +i := b k . These satisfy 
b k -\-1 Qfc+i "y^ (b k Rfc). 

Let fix = 1 — 4y and p 2 = 1 — 4^ with so chosen so small that 1 — 4^- > 2~i. This 
yields the conditions in (a) and (/3) and hence proves the desired estimate. □ 


3.3.2. Proof of Proposition 3.22. 

Proof of Proposition 3.22. The proof is standard (c.f. Theorem 6.9 of [PSU12]). 
We only give a brief sketch of it and indicate the main modifications. The idea is to 
apply the boundary Harnack inequality (Lemma 3.25) to the tangential derivatives 
d e w , e £ C;(e„). This gives the uniform (in h) C 1,a regularity of level sets {w; = h} 
with small h > 0. Then the uniform convergence of the level sets as h —> 0 allows 
to deduce the C 1,a regularity of the zero level set, which is the free boundary Tu,. 


We consider tangential derivatives v := d e w with e £ {x £ S n D B[\x n > r]\x"\}. 
Here p = C(n,p) max{eoj c»} is the same constant as in (29). Using the off-diagonal 
assumption (Al) and extending v to the whole ball B\ as described in Remark 3.9, 
we note that it solves 


did^djV = diF 1 in B\ \ A W: 
v = 0 on B[ (~l A w . 


By Lemma 4.3 in [KRS15], || dist(-, T^) - ^ ln(dist(-, r u ,)) _2 U*|| iP ( Bl ) < C(n,p)c*. 
Moreover, by (29), v satisfies the non-degeneracy conditions in Lemma 3.25 (after 
passing to the normalized function jf-). Thus, for eo and c* sufficiently small, 
the boundary Harnack inequality (Lemma 3.25) yields the existence of constants 
a £ (0,1), 9 = 9{n,p) and C = C(n,p) such that for all xq £ T w D B[/ 4 


(32) 


d e w d e w 

O) - a - 


d n w 


d n w 


< C^—{e n /2)\x - x 0 \ 
o n w 


x £ (xq + C' g (e n )) n Bi(x 0 )- 


We observe that for j £ {1,..., n — 1}, the estimate (32) can be improved to the 
following more quantitative control 


djW_ 

d n W 


(x) 


djW_ 

d n w 


(x 0 ) 


(33) 


< C'max{e 0 , c »}\x — a; 0 |“. 
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Indeed, given ej £ S n r\B[ with j £ l},we consider the vector e = ciej + 

C 2 e„, where the constants Ci,C 2 are chosen such that e £ {ir £ S n fi B[ | x n > rj\x''\}. 
We note that (for sufficiently small values of rj) it is possible to consider Ci > ^ and 
C{n,p) max{eo,c*} < C 2 < 2C{n,p) maxjeo, c*}. Inserting this into (32) thus leads 
to (33). 

We claim that (32) and (33) imply the desired C 1,a regularity of the free boundary 
and give the desired explicit bound for the Holder constant of g. Indeed, heuristi- 
cally this follows from the identity 


(34) 


dj9(x") = - 

o n w 


(a;",p(x"),0) 


and the fact that \b<\a Holder continuous up to Rigorously, we have 

to justify that the derivative of g exists and satisfies the claimed identity. In fact, 
since d n w > 0 in B[ \ A w , for h > 0 small, the level sets {u> = h} D B[ are given as 
graphs x n = gh(x "), where {gh}h>o is an increasing sequence in h which, by (32), 
is uniformly bounded in C' 1 (-Bj , ). Thus, gh converges uniformly to g as h goes to 
zero by Dini’s theorem. Next we notice that 


djg h (x") 



(x" ,gh(x"),0) 


Therefore, by (32), 


9jgh{x") - 



(x",g(x"),0) 


djw 

djw 


d n w 

(x",g h (x"), 0 ) 

(x",g{x"), 0 ) 


< C \ 9 h{x") — g(x")\ a . 


Since the sequence gh converges to g uniformly, we infer 


y 0 as h —^ 0. 

Therefore, by the fundamental theorem of calculus, the partial derivatives of g exist 
and they satisfy the identity (34). Moreover, by (32) we have djg £ C°'“(r ffl nBy 4 ). 
This finishes the proof of Proposition 3.22. □ 


dj 9 h(x") - ( 


(x",g(x"),0) 


Remark 3.26. Under the conditions of Proposition 3.22, the C 1,a norm of g is 
bounded by a constant, which only depends on n andp. Indeed, by Proposition 3.8 
and (33) (and by g( 0) = 0) we have that 

||Vg||cT.“( b" /2 ) < C{n,p) max{e 0 , c*}. 

Moreover, the proof of Lemma 3.25 and the proof of Proposition 3.22 show that 
in this case the Holder exponent a £ (0,1) only depends on n andp (if eo,c* are 
chosen sufficiently small in dependence ofn,p). 

3.4. Proof of Proposition 3.10. Last but not least, we come to the proof of 
Proposition 3.10. For convenience, we introduce the following notation 

d(x) := dist(x, T). 
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The proof of Proposition 3.10 is a consequence of weighted energy estimates 
(with weights which are naturally induced by the symbol of the operator) and local 
pointwise estimates. 


Proof. Step 1: Existence and uniqueness. We define a Hilbert space P in which we 
solve an appropriate weak formulation of our problem. Let u,v £ Cq°(R" +1 \ A). 
Consider the following inner product 

(■,•)«: Cg°(R n+1 \ A) x C“(R” +1 \ A) -a R, 

(■ u,v)u~ j Vu • \7vdx + J d{x)~ 2 uvdx for all u, v £ P. 
r+i 

Using this, we define our Hilbert space P as the closure of Co°(R n+1 \ A) with 
respect to ||m|| w := (u,u)u: 


H := C'o°(R"+ 1 \ A) H 

On P the bilinear form associated with the equation (19), 


B(u,v) = J 


a? 3 diudjV + K 


J d(x)~ 


uv, u,v £P, 


R ”! 1 8 "+! 

is bounded and coercive, thus by the Lax-Milgram theorem a unique solution exists 
in P. This proves the first part of Proposition 3.10. 


Step 2: Weighted energy estimate. Let g be the intrinsic metric determined by 
the symbol d(x) 2 |£| 2 in the cotangent or by the symbol d(a ;)~ 2 |£| 2 in the tangent 
space, i.e 

g x (v, u>) = d(x)~ 2 v ■ w, v, w £ T X (R " +1 \ T), x £ R " +1 \ T. 

The distance function associated with g is given by 

(35) dg(xi, X 2 ) = inf / d(x(s))~ 1 ds , 

c Jc 

where the infimum is taken over all rectifiable curves C in R n+1 \ T joining two 
points X\,X2 £ R n+1 \ T parametrized by arc length. We note that dg(-, •) is finite 
on each connected component of R ” +1 \ T. Let / : R n+1 \ T —> R be a Lipschitz 
function w.r.t. the intrinsic metric dg. It is not hard to see that / is Lipschitz w.r.t. 
dg if and only if it is locally Lipschitz w.r.t. the Euclidean distance and it satisfies 

|V/(s)| < Cd(x)~ 1 for a.e. x £ R ” +1 \T. 

We decompose the open set R n+1 \ T into Whitney cubes {Qj} with respect to 
the Euclidean distance in R ” +1 \T (c.f. Section 3.2.1 for the definition of a Whitney 
decomposition). Without loss of generality, we may assume that the decomposition 
is symmetric with respect to the x n +± variable. Indeed, the symmetry assumption 
can always be realized as Tc R" x {0}. 

We now make the following claim: 

Claim. For any k > 0, there exists s = s(n, k) > 0 sufficiently large, such that for 
any compact subset W C R n+1 \ T w have 

OO 

|| d -R e -s d s( x ’ w ) 


(36) 


< C'(n)dist(W,T) K . 
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We postpone the proof of the claim to step 4 and continue the proof of Propo¬ 
sition 3.10. 


Let <f> be a Lipschitz function w.r.t. the intrinsic distance dg and assume that 
\(/>(x)\ < re|lnd(a;)| for some k > 0. Consider u := e^u. Then (in a distributional 
sense) u is a solution of 

di(a l -’djU ) — {di4>)a 1 ^ djU — di(a l -’udj(j)) + {di4>){dj(ji)a^ u — Kd~ 2 u 

= d i (e^ i )-(3 i ^P i + eV 


We insert e~ sd ^ x,w ^u as a test function into the divergence form of (37), where 
s = s(n, k) is the constant from the claim and W C K. ra+1 \ T is a given set. Note 
that e - sd sO ’ w )u is a function in %. Indeed, using the definition of u, Holder’s 
inequality and the statement of the claim yields 

||d - le -„J B (x >W ) fi || La(RB+1) < ]^ d -K—l e ~sdg (*>W) u || i2(Rn+1) 

< sup ||d- 1 «|| La(<5) 

\ j J Qc{Q fc }fc 

< C(n) dist(W,r) _ ' t ||d” 1 u|| L 2 ( R n+i). 

Similarly is is possible to estimate V(e~ sds ^ x,w ^u). As a consequence of the ad¬ 
missibility of e~ sds ( x,w ' > u as a test function, we infer that if K = K(\\(f>\\ c o_,i,n) is 

dg 

large (where || • H^o.i denotes the Lipschitz semi-norm with respect to the metric 

dg 

dg), then there exists C > 0 depending on the cllipticity constants and dimension 
n such that: 


|| e -.« J ,(x > w) Vfi || La(RB+1) + K /2\\d- 1 e- sd ^ w ^u\\ L 2 ( ^ ) 
< C (|| e -^(xW) e ^i|| i2(Rn+1) + || e - s d 9 -(x, W ) e ^ 5 || i2(Rn+i) j _ 

This can be written in terms of u as 

|| e - sd9 (x,VE) e 0 Vu || i2(Rn+i) + || e - s d 9 -(x,iv) e ^-i u || i2(Rii+i) 


(38) < C (|| e -«*»(*.wO e ^|| i2(R „ +1) + ||e-^(-^) e ^ 5 || i2(Rn+1) 


)■ 


Step 3: Local pointwise estimate: Consider a given point x £ K" +1 \ T and let 
Q £ {Qj} be the Whitney cube which contains x. Let 


d := d(x). 


By the weak maximum principle (c.f. Theorems 8.15 and 8.17 in [GT01]) as well as 
by the properties (Wl) and (W2) from the definition of a Whitney decomposition 
(c.f. Section 3.2.1), there exists C = C(n,p) such that 


( 39 ) 


K*)l < c(d-^\\u\\ LHM . ) + d 1 - a J 1 ||*’ i || £ , (V<j) +d 2(1 -^ ) \\9\\L^ <3) ) , 
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where Mq is the union of all the Whitney cubes which touch Q. Recalling the 
assumptions on F l and g and using the properties (Wl), (W2), we have 

(40) II^IIlpM,) < 

(41) \\9\\lp/*{Mq) ~ IM 1 p "fflLp/ 2( r +1 )d V 1+ p ■ 

To estimate d~^ n+1 ^ 2 \\u\\L^(j^^), we seek to apply (38) with W = Mq and 
(fix) = —/tin d(x), k = (n + 1) ( -) + a. 

V 2 pJ 

The Lipschitz constant (w.r.t. dg ) of <f> is bounded by c n (a + 1). Thus, by taking 
K = C(er + 1) with C = C(n) large, we have (38). We notice that from the 
properties (Wl) and (W2) for Whitney cubes as well as (38), 

< d v+ 1 -^||d- K - 1 e-^^’^) M || i2(R „ +1) 

< d v+ 1 - 2 v i (||d-"e- sd 5(^W < 3 ) F i|| i2(Rn+1) 

+ \\d-« +1 e- sd ^’^g\\ L ^ + i ) ). 


Due to Holder’s inequality, property (W2) and due to the definition of k 


\\d-*F'\\ L . W) < \\d- a F i \\ LP(R n + l ) \\d- K+a \\ 2J, 

LP-J (Q') 

< || d °’T 1l || I/ p( R n + l), 

IM _k+1 .9||l 2 (Q') < Wd 1 ~ ct 9||lp/2(r"+i), 

for any Q' £ {Qj}- Thus, combining this estimate with (36) (applied with k = 0) 
leads to 


(43) 

\\d~ K e -sdg{x,MQ) pi 


L 2 (Rp+1) < ^ 


Wp-sdgixMo) 


sup ||d K F l || L 2 (q0 
dQj) Q'e{Qj} 


< C(n)\\d~ a F i \\ L p^ K n.+iy 


Analogously, we infer 

\\ d -^ e sM^) gh2(Rn+1) < C(n)||d 1 "^- CT 5 || iP/2(R „ +1) . 
Combining (39)-(41), (42) and (43), yields a constant Co = Co(n,p) such that 


\u(x)\ < Cod a+1 p (||d <t ^’ 1 ||lp(rp+ 1 ) + IM 1 p CT g|Lp/ 2 (Rp+i)) ■ 

Step 4 : Proof of the claim. Without loss of generality, we only prove the claim 
for W = M(Q). Let Mi := {Qj\ Qj touches M(Q)} be the set of all Whitney 
cubes which touch M(Q). In general, for k > 2 let A4 denote all Whitney cubes 
which touch at least one of the Whitney cubes in Mk-i- Hence A4 {Qj}- It is 
immediate from property (W3) for Whitney cubes that: 

(a) Let jjA/fc denote the number of the Whitney cubes contained in Mk- Define 
Wo = 0. Then Wk - Wk-i < (12) fc ("+ 1 ), k > 1. 
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Moreover, from (35), property (Wl) for Whitney cubes and an induction argument, 
it is not hard to see that: 

(/3) For any x € Qj € A 4, k > 2 we have dg(x,Af(Q)) > -A^==. 

Let r 0 be the minimal sidelength of Qj £ Af(Q). Note that by (Wl) and (W2) we 
have ro < dist(A/’ ( j, T). Combining (a) and (/ 3 ), we obtain 

OO 

'EW< rK e- ad ’< x ' U * ) \\ L - (Qi) 

i=i 

OO 

<r 0 -« + (4- 1 r 0 )-' ( |tA4 + X) E ^~ k ro)- K \\e- d ^^\\ L ^ Qj) 

k=2 Qj GA/" k\fsfk — i 

OO ( fe 1) 

< + 4 K 12" + Vo K + r(f K 12 fe(n+1) 

2 

< C'(n)r ( C K , 

if s = s(n, re) is chosen to be sufficiently large. 

Step 5: Continuous realization of the boundary data. In order to obtain the 
continuous attainment of the boundary data, we observe that any point x € A lies 
in an associated Whitney cube Q. Since the Whitney decomposition is considered 
with respect to T, the cube Q does not intersect T. As T = 9A, this implies that 
Q fl (R” x {0}) C A. Moreover, the same is true for all neighboring Whitney cubes 
Af(Q). Hence, we find a ball Bq which contains Q such that 

dia^djU — d(x)~ 2 u = diF 1 + g in Bq , 

u = 0 on Bq fl (R" x {0}). 

In particular, this implies that the boundary data of u are attained continuously 
in Q (c.f. [GT01], Corollary 8.28). As this holds for any point x G A, it proves the 
desired continuous attainment of the boundary data. □ 

Remark 3.27. In order to deduce the improved behavior close to A from Remark 
3.12, we argue by elliptic regularity and a scaling argument. In the case of a point 
xq £ B\ with d(X q) ~ 1 but dist(iroj A) < c n (where c n is a dimensional constant 
less than 1/2) we have the bound 

IMIwg.^Bi^Ozo)) < C{\\d~ cr F\\ L p( Rn +i' ) + \\d~ 1 ~ a g\\ L p/2( R n+i* ) ). 

By the Sobolev embedding theorem and the fact that u vanishes on A, this implies 
that 

|m(zo)| < C'(IM _ ' 7 ^||lp(r«+i) + |M 1-<T <?IIlp/ 2 (R"+ 1 )) dist(ico, A ) 1 . 

The general case follows by rescaling to the set-up with d(x o) ~ 1. 

4. Optimal Regularity and Homogeneity of Blow-Up Solutions 

In this section we return to the study of solutions of the thin obstacle problem 
(4) and present some consequences of the C 1,a regularity of the free boundary. Here 
our main result is an optimal regularity result for solutions of (4), c.f. Theorem 4, 
and a leading order asymptotic expansion of the solution at free boundary points, 
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c.f. Proposition 4.6. 

In the sequel, we assume that the metric a*- 3 satisfies the assumptions (Al), (A2), 
(A3), (A4). Moreover, the metric aU and the solution w are extended to Bf in 
the same way as in Remark 3.9. In this section we mainly work with the equation 
of djW in our C 1,a slit domain. For j £ { 1,... ,n}, djW satisfies (16) in B\ \ r„. 
For j = n + 1, we extend w to Bf by setting w(x', x n +i) = —w(x', —x n +i). Then 
d n+ \w £ C 0,ll/2 ~ s (Bi\fl w )nH 1 (Bi\fl w ) for any S > 0, as w = 0 on fl w . Moreover, 
it solves the divergence form equation 

did^djV = diF 1 in B\ \ fi TO , where F l = —{d n +\a?i)djW, 
v = 0 on fl w . 

Our main results rely on identifying the asymptotic behavior of solutions, c.f. 
Lemma 4.2 and Propositions 4.3 and 4.6, for which we invoke the C 1,Q regularity 
of the regular free boundary. This then implies the optimal local growth behavior 
(c.f. Corollaries 4.9-4.11) which in combination with our Carleman estimate from 
[KRS15] yields the optimal regularity of solutions. 

This section is organized as follows: We first deduce the behavior of general 
solutions of the Dirichlet problem in slit domains with C 1,a slits in Section 4.1. 
Then in Section 4.2 we apply these results to solutions of the thin obstacle problem, 
which yields the optimal local non-degeneracy and local growth estimates as well as 
a leading order asymptotic expansion in appropriate cones. In Section 4.3 the local 
results are then improved to global estimates and the optimal regularity result. 
Finally, in Section 4.4 we prove uniform lower bounds for a suitably normalized 
solution of the variable coefficient thin obstacle problem at regular free boundary 
points. 

4.1. Asymptotics and growth estimates. In this subsection we exploit the C 1,a 
regularity of the free boundary, in order to obtain a leading order asymptotic ex¬ 
pansion of solutions of elliptic equations with controlled inhomogeneities in C 1,a 
slit domains (Proposition 4.3). This expansion holds in appropriate cones. Using 
the asymptotic behavior, we then immediately obtain (local) growth estimates in 
the form of a Hopf principle (Corollary 4.4) and the optimal upper growth bounds 
(Corollary 4.5). Our approach relies on applying the boundary Harnack estimate 
to a given solution and a model function, which we construct in Lemma 4.2. 

In the whole of this section, we work under the following assumptions: 
Assumption 4.1. We assume that A and T are closed subsets o/R n x {0} with 
A = {x n < g(x")}, r = {x n = g(x")}, 

where g £ ^’“(R™ -1 ) for some a £ ( 0 , 1 ] and 3 ( 0 ) = 0 . 

As in the previous section, we use the abbreviation L for the operator L = 
dia^dj, and Lq for Lq = a^dfj. Again we use the abbreviations from Notation 
3.15 and set to := (2 4 v / n) _1 and c* := ||Va IJ 
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Under these assumptions we now refine the estimates from the barrier construc¬ 
tion from Proposition 3.16 in Section 3, in order to identify the leading order con¬ 
tribution in the solutions to our equations close to free boundary points. 

Lemma 4.2. Let A, T and g satisfy the conditions in Assumption 4-1 and let to, c* 
be as in Notation 3.15. There exists a function ho : B\ —> R such that 

(i) ho > 0 on B i \ A and ho = 0 on A, 

(ii) ho{x) > c n dist(:r, P ) 1 / 2 in B\ n {dist(a;, A) > £ 0 dist(z, T)}, 

(iii) Lho = 3 i+ 32 ; where g\ and 32 satisfy 

II dist(-,r)53i|| iP(Bi) < Cnc,, 

||dist(-,r)i- min <“’ 1 - Ii f i >32||L~(B 1 ) < C n (c* + ||V" 3 || c o. Q(Br) ) . 

Proof. As the proof follows along the lines of Proposition 3.16, we only indicate the 
main differences. 

Similarly as in our previous construction, we consider 

u°(x) := w 1 / 2 (x n ,x n+ i), 

and define ( x ) by composition with the change of coordinates T k as in Proposition 
3.16. Using these building blocks, we define h 0 in analogy to Proposition 3.16, i.e. 
for an appropriate partition of unity, r) k , we set: 

h 0 {x) ■.= '^2,g k vl(x). 

k 

By construction ho satisfies (i)-(ii). 

To show (iii), we first notice that Lh 0 = (dia l 3 )(djho) + L 0 h 0 := 31 + 32 - The bound 
for 31 follows immediately by noting that |Vh 0 (x)l < C n dist(a;, T) -1 ' 2 . In order to 
bound 32 , we compute Lo{ho) similarly as in step 1 of Proposition 3.16. Combining 
the fact that for each k, a *- 7 (x k )dijV®(x) = 0 in M ” +1 \ {x n +i = 0 , (x — x k ) • n k < 0 }, 
with the properties of the Whitney decomposition leads to 

y;(fl tJ (Xk)djjVk(x))Vk 0*0 = 0 in B 1 \ A. 

k 

We proceed with the estimate for the error contributions: Let v k and 17 denote the 
normals of A at x k £ F and xt £ T. Here x k ,xg £ P fl B\ are the projection points 
of two centers, x k and xe, of neighboring Whitney cubes, Q k and Qe. For these 
normals we infer \v k — U(\ < || g\\c°’ c ‘(B' 1 / ) r e ■ Thus, defining 7 := 1 — and 
following the computation in Proposition 3.16, we have that for all x £ Qn , 

(x k )dijr] k (x))v 0 k {x) < C(n ) (c*r7 5+7 + || Vg\\ c o.o,^)r ^' + “) , 

k 

5 Z^&^diVkdjvl < C(n ) (c*rf* +1 + \\V'g\\ C o^(B' 1 ')r^ +a ^ , 

k 

^2(a lJ {x) - a 13 (x k ))dij(r] k Vk(x)) < C(n)c*rf ' +7 . 

k 

Thus, the bound for 32 follows. □ 

We continue by deriving the boundary asymptotic behavior of solutions u to 
general divergence form equations in slit domains with C 1,a slit. The main tool is 
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our boundary Harnack inequality (Lemma 3.25). 


Proposition 4.3. Let A,F and g be as in Assumption 4-1. Suppose that u £ 
C(Bi) D is a solution of 


(44) dia^djU = diF 1 + ji in Hi \ A, u = 0 on A. 


Further assume that a lJ , F z , gi and u satisfy the assumptions from Proposition 3.17. 
If c*, S and || g\\c°’ a (B^) are sufficiently small depending on n, p and a, then there 
exists a constant (3 £ (0,min{a, 1 — a C 0,/3 (r n S 4 / 2 ) function b : T —> K., 

b > 0, such that for all xq £ T fl B[j 4 and x £ B 1 / 4 ( 20 ) 


(45) 


u(x) - b(xo)w 1 / 2 


V Xo • ( X - Xq) X n+ 1 \ 

,K • A{X 0 )V X0 Y/^ {a n+l,n+l {xo) )l/2) 

< C ^c* + ||V ,, 5 || C r 0l a ( B »\x — xo \ 2+ ^. 


Here u Xo is the outer unit normal of A at xo £ F , g denotes the parametrization of 
r and A(x) := (a lJ (*)). 


Proof. We seek to apply the boundary Harnack inequality (Lemma 3.25) to the 
pair it, ho (with <5o = min{a, 1 — ^^±1}), where h 0 is the function from Lemma 4.2. 
By Lemma 4.2, ho satisfies the conditions in Lemma 3.25. Moreover, we note that 
by Proposition 3.17 

(46) u(x) > c n dist(x, T)5+f in B 1 / 2 Fl {dist(x, A) > £ 0 dist(x, F)}, 
and 


(47) 


u(x) > — dist(x, T) 2+6 in B 1 / 2 


Here e = a + £ — and a is the constant from Proposition 3.17. Hence, by the 
boundary Harnack inequality we have that 


(48) 


sq- 


u(e n /2 ) 


< 


u(x) _! u(e n /2 ) 
h 0 {x) ~ h ° h 0 (e n /2) ’ 


^o(^n/2) rio\x) 

for x £ H!/ 2 n{dist(a;, A) > £0 dist(x, F)} and for some So = so(n,p, || S7"g\\co,c 
Moreover, ^( 3 .) is Holder continuous in Cg(rnH(y 4 ) up to FnH^ 4 . More precisely, 
for each xo £ rni? 4 ^ 4 , the limit of as x —> xo, x £ xo +Cg(e n ), exists. Denoting 
it by b(xo), gives 


(49) 


u(x) 


h 0 (x) 


- b(x 0 ) 


x /-t u ip n /2) | ,3 

° S ’ 


x £ ( Xq + Cg(e n )) Fl Bi/ 2 , 


where C = C(n,p ), for some 0 £ (0,1). By (46), Lemma 4.2 (iii) for ho and by the 
fact that g £ C 1, “(H"), we have u(e n /2)/ho(e n /2) > c n > 0. Combining this with 
the lower bound in (48) implies that b(xo) > 0. 

We now extend this to the whole neighborhood Bi/ 4 ( 2 : 0 ) instead of only considering 
the cones xq + Cg(e n ): Using (47) on all scales and arguing as in the proof of 
Lemma 3.25, we obtain for x £ B l / i {x 0 ) 


u(x) - b{x 0 )ho{x) > -C\x - l0 | mi ' , {|+^+« j 
u(x) - b(xo)ho{x) < C\x - a;o| min {l+ ? -3+«. 
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Consequently, for (3 = min{/3, e}, 

(50) |u(x) ~ b(x 0 )h 0 (x)\ < C\x — Xq\^ + i3 for x £ -Bi/ 4 (xo)- 

Next we show that for each x 0 € T fl By 4 , 

(51) 

\h 0 (x)-w 1/2 (T Xo x)\ < C (c* + \\^"g\\ c «(B'; /2 )) l*-a ; o |’ +min{a ’ 1 x £ B 1/2 , 

where T Xo is defined as the following affine transformation in R" +1 : T Xo (x) = 
R Xo B~^(x — xq), where = A(x o) and R Xo £ SO(n + 1) is such that 

Rx 0 Bl. 0 v X0 = c 1 >xo e n , R xo B t Xo e n+1 = c 2 , X0 e n+1 . A direct computation leads to 

(ci,x 0 ) = (RxoB Xq V X q ) ' {Rxa^x 0 v xo) = v x 0 ' A{xq)v Xo i 
(c 2 ,x 0 ) 2 = {R X0 B t X0 e n+ 1 ) • {R Xo Bl o e n+ 1 ) = a n+ 1 ’ n+ 1 (x 0 ). 

In order to show (51), we note that by a similar estimate as (22) with s = 0, we 
deduce 

(52) 

\wi/ 2 {T k x) - w 1/2 {T Xo x)\ < C (c* + ||V" S || c «, (b „ 3) ) \x - x 0 |5 +mm{a ’ 1 "^ } , 
for all x £ Q £ J\f(Qk)- Instead of Lemma 3.14, we use 

K - "xo\ < C || V"5|| c o,<, (s // ) \x k ~ *o| a , 

in this context (we recall that is the normal of T at Xk, where Xk £ F realizes 
the distance of the center Xk of Qfc to T). 

Thus, by using the properties (W1)-(W3) of the Whitney decomposition, we have 

\ho(x) - w 1 / 2 {T Xq x)\ < ^2 Vk\wi/ 2 (Tkx) - w 1/2 (T Xq x)\ 
k 

< C (c* + \\V'g\\ ca{ B^ \x x 0 |i +min{ “’ 1 _: 2 J 1} . 
Writing T Xo out explicitly, we therefore obtain 

v xo ■ (x-x 0 ) x n+1 \ 

K^(®oK) 1/2 ’ (a"+^(x 0 )y/2) 

< c (c* + ||V w 5 || c< , (B „ a) ) \x - solH-W- 2 ?}. 

Using the above asymptotics of h 0 in (50), we arrive at the desired result. □ 

The above asymptotic estimate leads to two immediate corollaries: On the one 
hand, it yields a Hopf principle and on the other hand, if combined with an interior 
elliptic estimate, it results in an upper growth bound. 

Corollary 4.4. Assume that the conditions of Proposition f.3 are satisfied and 
that 0 £ T. Then, if c*,<5 and || ’V"g\\c°’ a (B") are sufficiently small depending 
n,p and a, there is a dimensional constant c n > 0 such that 

u{x) > c n dist(x,T )2 

for x £ By 2 x (—1/2,1/2) fl (x| dist(x, A) > £ 0 dist(x, T)}. 


h 0 (x) - w 1/2 


on 
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Corollary 4.5. Assume that the conditions of Proposition 4-3 are satisfied. If c*, <5 
and || Hco, a (s'') are sufficiently small depending on n,p and a, then there exists 
a constant C = C(n,p) > 0 such that for all x £ B 1 / 2 

|u(x)| < C'M 0 dist(x,r)^, 

where Mq := sup Bl |u|. 

4.2. Local 3/2-growth estimate, asymptotics and blow-up uniqueness. In 

this section we transfer the results from the previous section into the setting of the 
thin obstacle problem. Hence, this yields leading order asymptotics in appropriate 
cones, a Hopf principle and an optimal local growth estimate. Last but not least, 
we explain how the asymptotics and the growth estimates can be used to derive 
the uniqueness of the 3/2 -homogeneous blow-ups. 

We begin by transferring the general results of Section 4.1 to the setting of 
solutions of the thin obstacle problem: 

Proposition 4.6 (First order asymptotics). Let w be a solution of (f) satisfying 
(AO), (Al), (A2), (A3), (Af). Assume that for some small positive constants eo 
and c* 

(i) ||w-wj 3/ 2 || c , 1(b + ) < eo, 

(ii) l|Vab|| iP(B + ) < c *. 

Then if e 0 and c* are chosen sufficiently small depending on n,p, there exists a 
function b e : n B 1 / 2 —> R, b e £ C°’ a for some a £ (0,1 — (where a can be 

chosen the same as the Holder regularity ofT w in Proposition 3.22), such that for all 
Xq £ T w nB 1/2 , we have the following asymptotic expansion for d e w, e £ S ,n_ 1 x{0}: 


(53) 


%n +1 


dMx) - b ’ Mwt/ 2 to ))'/ 2 


< C(n,p) max{eo,c*}dist(x,r w )5+ a ) x £ ^^(xo). 
Here v XQ is the outer unit normal of A w at xq and ^4(xo) = (a y (xo)). 

Before coming to the proof of this result, we make the following observations: 


Remark 4.7. (i) Similarly as in (53), it is possible to obtain an asymptotic 

expansion for d n +iw: There exists a function b n + 1 £ C ,0,a (r iu ), such that 
for xq £ Tu, fl B(i 2 and x £ B i/ 4 (xo) 


d n+1 w(x) - b n+1 (x 0 )w 1 / 2 


(x - x 0 ) • v Xo x n+ i \ 

K 0 • A{x 0 )u xo )H2' {a n + l,n+l {xo)) l/2 ) 
< C(n,p) max{e 0 , c*} dist(x, ru,)^ -1 "®, 


where w 1 / 2 (x) = — Im(x n + ix n +f) 1 ' 2 . 

(ii) Later, we will see that b e and b n +\ satisfy the following relation (c.f. (57) 
and (58) in Corollary f.8): There exists a C 0,a function a : —> K, 

a > 0 , such that 


b e {x 0 ) 


3 q(xp)(e • Vx 0 ) 

2 (y x o • A(xo)v Xo ) 1 l 2 1 


b n+ i(x 0 ) 


3_ a(xp) _ 

2 ( a n+r>n+I( xq)) 1 / 2 ’ 
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Proof of Proposition f.6. The proof of Proposition 4.6 follows immediately from 
Proposition 4.3. More precisely, we first note that, by Proposition 3.22 if c*,eo are 
sufficiently small depending on n,p, then (up to a rotation) 

r™ n b 1 / 2 = {x n = g{x")} n b 1/2 , 

where g £ C 1,a , ||si|<ji,a < max{c*, eo} and a only depends on n and p. Let v = d e ui 
and recall from the proof of Proposition 3.8 that it satisfies the following divergence 
form equation in B± \ A w 

did^djV = diF ", F l = —(d e a i: ’)djW, 

with || dist(-,r w ) _ 5 ln(dist(-, r u ,)) _ 2 F*|| i p( Bl ) < C(n,p)c » (c.f. Lemma 4.3 in 
[KRS15]). Moreover, d e w satisfies the assumptions of Proposition 4.3 if e is in 
a sufficiently large cone of directions, i.e. e £ C'(e„) with rj > 0 sufficiently small. 
Thus applying Proposition 4.3 we obtain the desired result. Note that the expo¬ 
nent a £ (0,1 — can be chosen to have the same value as the Holder regularity 
exponent of in Proposition 3.22, since both of these come from the boundary 
Harnack inequality (Lemma 3.25). The result for general directions e follows from 
Remark 4.7 (ii). □ 


As a corollary of the asymptotic expansion for Vw, it is now possible to obtain 
the asymptotic expansion for w via integration along appropriate contours: 

Corollary 4.8. Under the assumptions of Proposition f.6, there exists a Holder 
continuous function a : T^, —> R, a > 0, such that for xo € r w n an d 
x £ B 1 / a {xq) with x — Xq £ span{^ 0 , e n+ \} 


(54) 


w(x) - a(x 0 )w 3 / 2 


•^n+1 


(x - x 0 ) • v Xo 

v X0 -A{x q )v X0 ) 1/2 ’ ( a Ti + 1 >™+ 1 ( x 0 )) 1 / 2 / 

< C(n,p) max{e 0 ,c*}dist(a;,r !(J ) 2 +a . 


Here v Xo is the outer unit normal of A w at Xo and A(xq) = (a^(x o)). 
In particular, w(xo + rx)/r 3//2 has a unique blow-up limit as r —> 0. 


Proof. The asymptotic expansion of w follows from the asymptotics for dv XQ w, 
d n +\w and an integration along an appropriate path in the plane spanned by 
u Xo , e n +i. Given xq £ T^,, we write 

x = x 0 + t\v Xo + t 2 e n+ 1 for all x £ span{i/ Xo , e n+ i}. 

In particular, t\ = (x — xo) ■ v Xo and t 2 = (x — Xq ) • e„+i = x n +i- We consider the 
restriction, w, of w onto the two-dimensional plane spanned by v Xo and e n +i, i.e. 

w(ti,t 2 ) :=w{x 0 +hu Xo +t 2 e n+ i). 

Then, 

d tl w(ti,t 2 ) = d„ XQ w(x o + t\v Xo +t 2 e n+ i), 
d t2 w(ti,t 2 ) = d en+1 w(x o + t\v XQ +t 2 e n+ i). 

From the asymptotics derived in Proposition 4.6, we obtain 

d tl w(ti,t 2 ) = b Vx w 1/2 {cfH 1 ,c 2 h 2 )+g 1 (t 1 ,t 2 ), 
dt 2 w(ti,t 2 ) = 6 „ + iM} 1 / 2 (c)" 1 ti,c^ 1 t 2 ) + g 2 (ti,t 2 ). 


( 55 ) 
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Here 


\gi(ti,t 2 )\, \g 2 {h, t 2 )\ < C{n,p) max{e 0 , c*}(|fi| 1/2+ “ + \t 2 \ 1/2+a ) 


and b v = b v (xq), b n+ 1 = 6 n +i(a:o)- The constants c i, c 2 are (for fixed x 3 ) 
defined as 

ci = {v xo ■ A{x 0 )v Xo )*, c 2 = {a n+1 ’ n+1 (x 0 ))^ ■ 

Since §(w 1 / 2 , tDi/ 2 ) = Va) 3 y 2 , we can rewrite (55) as 

(56) 


25,, ci . . 

d tl w(ti,t 2 ) = —x 2 —5 tl 103 / 2(03 h,c 2 1 t 2 ) +gi{h,t 2 ), 


Since 


we have 


Here 


d t2 w(t!,t 2 ) = 2K + lC2 d t2 w 3/2 ( Cl 1 t 1 ,c 2 H 2 ) + g 2 (ti,t 2 ). 


j d tl wtd t2 (f> = J d t2 wd tl 4> for any <fi £ CZ°(R 2 ), 


R 2 


R 2 


^ 3 ° 1 ~ 26 "+ lC2 ^ J w 3/2 (c 1 1 t 1 ,c 2 1 t 2 )d? lt2 cl) = J g(t 1 ,t 2 )d? lt2 (j). 

R 2 R 2 


rti rt 2 

g{ti,t 2 ) = gx (*i,t 2 ) - g 2 (ti, t 2 ) = - / gi(s,t 2 )ds+ g 2 (t lt s)ds. 

Jo Jo 

A direct computation then gives 

\g{ti,t 2 )\ < C'(n,p)max{e 0 ,c»}(|ti| 3/2+ “ + |ti||t 2 | 1/2+Q 

+ \t 2 \\h\ 1/2+a + \t 2 \ 3/2+a ), 

which is of higher (than 3/2) order in (fi,f 2 ). Thus necessarily we have 
(57) 6 ^ 0 ci = b n+1 c 2 . 

Now let £(s) be the path from (0,0) to (ii,t 2 ): £(s) = (fis,f 2 s), s £ [0,1], and let 
£(s) = (ci 1 t 1 s,c 2 1 t 2 s). Then 

< C{n,p) max{e 0 , c*}|/(s)| 1 / 2+a . 


d w(£( S ))~^±^lA. w 3 /2 (£ (s » 


ds 3 ds 

Integrating from 0 to 1 leads to 

~ /, 1 \ ‘^bn-\-lC 2 / _ 1 , — 1 , \ 

W\tl,t 2 ) --- U>3/ 2 {C 1 h,C 2 t 2 ) 


Finally defining 
(58) 


i(zo) := 


<C'(n,p)max{e 0 ,c*}(|t 1 | 3 / 2+ “ + |t 2 | 3 / 2+ “). 
26 n+ i(xo)c 2 


and recalling the explicit expression of Ci and c 2 yields the asymptotics (54) for w. 
The fact that a > 0 follows from Corollary 4.9. □ 


From Proposition 4.6 and Remark 4.7 (i), we immediately obtain the following 
lower and upper bounds on the growth rate of d e w away from the regular free 
boundary: 
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Corollary 4.9 (Lower bound). Under the assumptions of Proposition f.6 there 
exist a constant Co = Co(n,p) > 0 and a radius ro = ro(n,p) such that 

d e w(x) > co dist(x, r^) 5 , x G {x G B^ 2 \ dist(x, A^,) > lo dist(a;, T^)}, 

for any e G {x G S n ~ l x {0}| x n > ||a;"|}. 

Remark 4.10. A uniform lower bound may not hold for d e w in a full neighborhood 
of the origin. However, after subtracting an “error term” or under higher regularity 
assumptions on the metric a lJ , it is possible to obtain a uniform lower bound in a 
full neighborhood of the origin (c.f. Section 4-4)- 

Corollary 4.11 (Upper bound). Under the assumptions of Proposition 4-6 there 
exist a constant C = C(n,p ) and a radius ro = ro(n,p) G (0,1) such that for each 
xq £ r w fl Bi /2 

(59) sup |Vw| < Cr 1 / 2 , 0 < r < ro- 

B r (x o) 

Remark 4.12. The above upper bound holds for solutions which are sufficiently 
close to W 3/2 (c.f. assumption (i) of Proposition 4-6). In the next section we 
will remove this assumption and show a uniform upper bound for solutions w with 
IMIl 2 (b +) = 1 (c.f. (60);. 

4.3. Optimal regularity of the solution. In this section we exploit the com¬ 
parison arguments from the previous section in combination with the Carleman 
estimates from [KRS15], to obtain the optimal regularity of solutions to the thin 
obstacle problem with coefficients aU G W 1,p with p 6 (n + l,oo]. In comparison to 
the results of Proposition 4.8 in [ vRS15] our main improvement here is a uniform, 
optimal upper growth bound without losses (c.f. (60)). This then allows us to 
remove the logarithmic losses from the regularity estimates from Proposition 4.8 in 
[KRS15] (which were exemplified in the non-uniform constant C( 7 ) > 0 and the 
only C 1,1 ' 2-6 regularity of solutions to the variable coefficient problem in Proposi¬ 
tion 4.8 in [KRS15]) and to obtain the following optimal regularity estimates: 

Theorem 4. Let w be a solution of (4) in B(/ which satisfies the normalization 
condition (AO) and let a 12 : B+ be a W 1,p tensor field with p G 

(n+l,oo] satisfying (Al), (A2), (A3) and (A4). Then the following statements 
hold: 

(i) There exist a constant G = C(n,p ) and a radius Ro = Ro(n,p, ||VaU 
such that 

(60) sup \w\ < Cr 3/2 for any r G (0, Rq) and Xq G D B\ . 

B r (xo) 2 

(ii) If p G (n + 1,2(n + 1)], it holds that w G C 1 , ' y (By 2 ) with 7=1 — 2±I. 

Moreover, there exists a constant C = C(n,p,\\a l I\\ wlp ^ B+ ^) (which re¬ 
mains bounded as 7 1 / 2 ) such that 

IVu^o; 1 ) — Vu>(:r 2 )| < C\x x — x 2 | 7 for all x x ,x 2 G By 2 . 

(iii) If p G [2(n + l),oo], then w G C' 1 , 5 (.B+ ). More precisely, there exists a 
constant C = C(n,p, \\a v : ||yyi,p(B+)) such that 

_ rr2 1 !/ 2 f^ n ll rr - 1 T 2 c n + 
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Proof. We will only show the growth estimate from part (i). With (60) at hand, 
arguing similarly as for Proposition 4.8 in [KRS15], we then obtain the correspond¬ 
ing optimal regularity results (ii) and (iii). 

By Lemma 4.1 in [KRS15] and by the gap of the vanishing order (i.e. either 
k Xo = 3/2 or k Xo > 2), it suffices to show (60) for xq £ r 3 / 2 («;) D B[j 2 . For sim¬ 
plicity we assume that 0 £ r 3 / 2 (ii') and show the growth estimate (60) at xo = 0 . 

Step 1: Suppose that w is a solution to the thin obstacle problem which satisfies 
the assumptions (AO), (Al), (A2), (A4) with a metric a y £ W 1 ’ p (Bf) for some 
P > Po, Po £ {n + l,oo). Suppose that 0 £ r 3 / 2 ('wi). We will show that for any 
6 > 0, there exists a constant e = e(n,p,po,S) > 0, such that if 

(i) ||a«-^|| Loo(Bi+) <e, 

(ii) ||® ■ Vw- §H| i3(A + 2ii) < e, 

then there exists a nontrivial 3/2-homogeneous global solution wo such that ||w — 

w o\\c 1 (B+ i ) - 3- 

Suppose that this were wrong, then there existed a parameter So > 0, a sequence 
€k —t 0 and a sequence, Wk, of solutions to the thin obstacle problem 

did^djWk = 0 in B+, 

Wk > 0, ~d n+ iw k > 0, w k d n+ iw k = 0 on B[, 

which satisfy the assumptions (A0), (Al), (A 2 ), (A3), (A4), such that 

(i) a/? £ W 1 ’ P (B +) for some p > p 0 with p 0 as above, i.e. po £ (n + 1, oo] and 
Po being independent of k, 

(ii) \Wk - Sij \\L°°(B+) < 

(iii) \\x ■ Vw k - |iWfc|li, 2 ( A + j < £fc, 

but Urcfe — wo\\ C i( B + ) > So for any nontrivial 3/2-homogeneous global solution wo- 

By the fundamental theorem of calculus, there exists an absolute constant C > 0 
with 


IKII £ 2 (A+ gjl ) 


<C\\w k \\ L 2 (A + ) + C|||a;| 3/2 a;-V(|x| 3/2 w fe )|| i2(j4 + } 

'-1/2,7/8' ^ l/2,l' 

- C\\ w k\\ L a(A+ /a 7/g ) + Cek\\wk\\ L 2( B +) 

— C\\wk\\ L 2(Bf /s ) + C e fc|| u; fc|| L 2 (A+ /g + Cek\\Wk\\ L 2( B +^y 


Hence, if £k > 0 is sufficiently small, then 


ll w fcllz, 2 (A+ gl ) - 2C'||tnfc|| L 2( B + g ), 

which implies that 


(61) 1 - 11^-11^2(5+) < (2C + l)\\wk\\ L 2 ( B +^y 

By the interior H 1 estimate and by Proposition 4.8 in [KRS15], Wk is uniformly 
bounded in B[ 1 (B^ g ) and C 1 ’' l0 (By s ) for some 70 = 7 o(n,po) > 0- Thus, up to a 
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subsequence 

w k — y w weakly in B: 1 (By 6 ), 
w k w strongly in L 2 (B+ /g ), 
w k ^ w in C\ oc {B+! g ), 

where w (by the contradiction assumption (i)) solves the constant coefhcient thin 
obstacle problem 

Aw = 0 in By s , 
w > 0, — dn+iw > 0, wd n +iw = 0 on By 8 . 

Furthermore, by the contradiction assumption (ii), w is homogeneous of degree 3/2 
in By g . By analyticity, w is a 3/2-homogeneous global solution in R” +1 . Moreover, 
by (61) and by the strong convergence in L 2 (By s ), 

IMk 2 (B 7/8 ) > (2C + 1) \ 

thus w is nontrivial. Therefore, we have found a nontrivial 3/2-homogeneous global 
solution w such that (up to a subsequence) ||uy — w\\qi^ b + ) —> 0 as k —> oo. This 
is a contradiction. 


Step 2: We show that there exists e = e(n,p), such that if for some r £ (0,i?o) 
with Rq = Ro(n,p, ||Vo* : '|| I(1 , (B + ) ) 


3 

(62) Ik' Vw- ^w\\ L 2(A+ /2r ) < c IMIl 2 (b+)> 

then there exist constants po = po{n,p) G (0,1) and C\ = Ci(n,p) > 1 such that 

H4I < \\w\\l*(b+) \ x \* for x e B tor- 

For this we first notice that (62) can be rewritten in terms of L 2 normalized 
rescalings 


(63) 


w r {x) := 


j(rx) 


r " 2 1 |M 


L 2 (B+) 


to yield 

3 

(64) Ik - V-uv — -uy|| r2M + 1 < e for some r £ (0, i? 0 )- 

2 '•1/2,1-' 


Here w r is a solution to the thin obstacle problem with coefficients a l J{x) := (rx), 

which satisfy ||V4 J '|| LP(B + ) = r 1 ”^ \\^a ij \\ LnB +y 

Secondly, by Proposition 3.22 in Section 3, there exist constants eo = eo (n,p) and 
c» = c*(n,p), which are less than some universal constant, such that if Huy — 
w 3/2\\c 1 (b + ) < eo with W 3 / 2 {x) = c n Ke(x n + ix n + 1) 3 ^ 2 (possibly after a rota¬ 
tion) and ||Va 1 - 7 || iP ( B + j < c„, then the free boundary of uy is a C 1,a graph in 

By T Moreover, its C 1,a Holder norm is uniformly bounded, depending only on 
n,p. The Holder exponent, a, only depends on n,p (c.f. Remark 3.26). Thus, by 
Corollary 4.11, there exist constants po and Ci depending on n,p, such that 

|uy(x)| < Cik| 3/2 for any x £ B+ 0 . 
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Scaling back then yields 

(65) Hz)| < C-ir-ir-^ |MIl 2 (b+) M* for any x G B+ or . 

To complete the proof of step 2, we apply step 1 to w r with 6 = e o- This deter¬ 
mines the parameter e = e(n,p) (if Rq = Ro(n,p, || Va 1 - 7 || iP ( B +)) is chosen such that 

ll Va r' 7 llL!>(B 1 +) < c*(n,p) for any r G (0,R o ])- 

Step 3: Fix e > 0 as in step 2. Let r\ G (0, i?o] be the largest radius such that 
(62) holds (we remark that the existence of an n > 0 such that (62) is satisfied, 
follows from the fact that the L 2 rescaling has a 3/2 homogeneous blow-up along a 
certain subsequence, c.f. Proposition 4.5 in [KRS15]). Then for any r G [ri,i? 0 ], 

(66) \\x • Vic - 3 /2w\\ L 2 (Ar/2 r) > e||w ;|| L2(A + ^ 

In this regime we discuss two cases: 

3 | Ti + 1 | e 

• Case 1 : |MIl 2 (b+) < r i 2 2 - hi this case 

b^i" 4 " IMIb2 ( b+) < rf. 

Thus, by (65) in step 2, there exists C = C(n,p , ||Vo® j '||i,p) such that 
|w(z)| < Cr?\x\ 3/2 for any x G B+ ri . 

Recalling the parameter dependence of po, we obtain 

(67) \w{x)\ < C\x\%+i for any x G A+ oruri . 

For any r G (ri,i? 0 ), we then use Corollary 3.1 in [KRS15] (note that this 
estimate has a logarithmic loss so that it is necessary to use the slightly 
improved estimate (67) in the application of the Carleman inequality; al¬ 
ternatively it would have been possible to use the Carleman estimate from 
Lemma 3.2 in [KRS15], for instance in the version of Remark 10 in [KRS15] 
which does not need the e-improvement). Thus, we obtain 

IMIl 2 (a+ /2 ) ^ Cri +r ^ with C = C(n,p, llVa^H) 

for all r G (ri, Rq). 

3 _|_ ra + 1 _|_ _e 

• Case 2: |H|.l 2 (.b+ ) > r/ 2 2 . In this case statement (i) from Lemma 

4.2 in [KRS15] (with the two radii n and Rq) would either imply that 
r\ is bounded from below, i.e. r\ > C , for some constant C > 0 which 
only depends on n,p, | V11 f p . or it would imply a contradiction to the 
maximal choice of r\. 

In summary, we either have a uniform lower bound for rr (depending only on 
n,p, || Va®- 7 || iP ) or we obtain that the growth estimate (60) holds for all r G (0, Rq). 

Combining steps 1-3, we have thus obtained a radius Ra = Ro(n,p , ||Va I:, || iP ) 
and a constant C = C(n,p) such that sup Br |u;| < Cr 3 / 2 for any r G (0, Ro). This, 
together with Proposition 4.8 in [KRS15], completes the proof of Theorem 4. □ 
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4.4. Improved lower bounds. In this section we improve our lower bounds from 
Corollary 4.9 by making them uniform: While in general a uniform version of 
Corollary 4.9 does not hold for the full function w, after a suitable splitting, it is 
possible to show that the leading order term satisfies a uniform lower bound. More 
precisely, let : B+ —*■ be a uniformly elliptic W 1,p , p £ (n + 1, oo], 

metric, let f £ L p (Bj) and let w be a solution of 

dia^d^w = f in Bf , 

(es) j J y 

d n +iw <0, w > 0, wdn+iw = 0 in B [, 

satisfying the normalization conditions (A0)-(A4) from Section 2 and conditions 
(i), (ii) from Proposition 4.6. 

As in Remark 3.9, we extend w, the metric a u as well as / from to B\ \ A w 
and further to K n+1 . We now split w into two components 

w = u + u, 

where u solves 

(69) a^dijU — dist(x, T w )~ 2 u = / — (dia^)djW in M” +1 \ A w , u = 0 on A w , 
and the function u solves 

(70) a^dijU = — dist(x, r^,) -2 ?/ in K" +1 \ A„j, u = 0on A w . 

As explained in the passage following Remark 3.9 in Section 3.2 the intuition is that 
u is a “controlled error” and that u captures the essential behavior of w. Moreover, 
as we will see later, u is C 2 ’ 1 p~ (K) regular for any K (e R n+1 \ r„,. 

Lemma 4.13 (Positivity). Let f £ L P (B+), cW £ W 1,p with p £ (2 (n + 1), oo], 
and u,w be as above. Then we have that u £ C 2 ' 1 p~ (K) D C 1 ’^ (B 1 / 2 ) for any 
K <s Bi\ r,,,. Moreover , for e £ C^ie-n), d e u satisfies the lower bound 

d e u(x) > cdist(x, A w ) dist(x, Tu,) -3 , c > 0. 


Remark 4.14. If f = 0, the statement of Lemma f.13 can be improved to hold for 
metrics a u £ W 1,p with p £ (n + 1, oo]. 


Proof. Let u be defined as in (69). By the refined estimate from Remark 3.12, we 
infer that 

(71) ... ' 1 -~ 


|u(x)| < c* dist(x, A^) dist(x, r m ) 


Moreover, combined with the C ’ p regularity of u away from r„,, we obtain 

the up to r^, regularity of u, i.e. u £ C ' p (BP and IImII j 1 _n±i < Cc*. 

C ’ v (B{) 

Thus, applying the classical elliptic estimates to u away from we deduce that 
u £ C ’ p (K) for any K <s B i \ T^. Recalling the decay estimate of w (c.f. 
Theorem 4) and combining it with the estimate (71) for u, yields that u = w — u 
satisfies 


(72) 


|u(x)| < Co dist(x, r w ) 2 . 


This together with the W 3,p interior estimate gives that u £ C X ’^(B\). Next we 
show the lower bound of d e u. This follows from an application of the comparison 
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principle from Proposition 3.17. First we note that at B\ D {|x n +i| = £o}, where 
to = -j= > 0, d e u is non-degenerate in the sense that 

(73) d e u(x) > c„ dist(x, T^)^, eeC' v (e n ). 

Indeed, for x G -Bin{|x ra+ i| = £o} we have dist(x, T^) ~ dist(x, A^), where by (71) 
the function d e u satisfies 

(74) \d e u(x)\ < c»dist(a;,r ii; ) 1__ p _ . 

On the other side, the asymptotics from Proposition 4.6 yield that d e w(x) > 
cdist(x, T^) 1 / 2 for x G B\ n {|x„+i| = to}- As p > 2(n + 1), dist(x, T^,) 1 / 2 domi¬ 
nates and thus d e u inherits the lower bound of d e w on x G B\ D {|x„+i| = £q}- 
Next we consider the equation of d e u (in non-divergence form): 

a^dijdeU = —( d e a l i)dijU — <9 e (dist(x, r^) -2 !!) 

=:H + diG\ 


The functions H, G l satisfy 

dist (x,T w )?H e L p {Bi) 7 
dist(x,r^) 11 ^ 1 G l G L°°(Si). 

Here the first estimate follows from the W 1,p regularity of a lJ and the pointwise 
interior estimates for dijW, the second estimate is a direct consequence of the previ¬ 
ously derived bound (71) for u. Furthermore, in B i we have d e u = d e w—d e u > —c *. 
Thus, the assumptions of Proposition 3.17 are satisfied (note that Proposition 3.17 
also holds for non-divergence form equation) and we obtain the desired lower bound 
for d e u in the region {x G 5i| dist(x,T™) ~ dist(x, A^,)}. 

In the end we remark that this lower bound holds in the whole ball. Indeed, this 
follows once more from a comparison argument: Let h = c s h s + ho, where h s is 
the barrier function from Proposition 3.16 and ho is the barrier from Lemma 4.2. 
Choosing s = s(a,n,p) appropriately yields that h satisfies the conditions (i) and 
(ii) of Proposition 3.17. Moreover, by using an up to A^ estimate for q(x) as well as 
the fact that ho dominates h s (for an appropriately chosen constant c s ), we obtain 

(75) h(x) > cdist(x, A u ,)dist(x,r uJ ) _ 5. 

Relying on the test function 

u(x) := d e u(x) + \x' — ( Xo)'\ 2 — 2 ~ s h,(x) — 2?rx 2 +1 , 

and introducing a similar splitting as in Proposition 3.17 then yields 

— n + 1 

d e u(x) > ch(x) — c* dist (a;, A^,) dist(x, r w ) p . 

The second term on the RHS originates from the estimate of the auxiliary function 
u\ (from a similar splitting as in the proof of Proposition 3.17), where we have 
exploited the Lipschitz regularity of G l away from r w (indeed, without this obser¬ 
vation, the error estimate would not suffice to absorb the error contribution, c.f. 
Remark 4.15 below). Recalling the lower bound for h (c.f. (75)) implies the desired 
lower bound for d e u. □ 
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Remark 4.15. We note that without the Lipschitz regularity of G l , we would only 
have had error estimates of the form 


d e u(x ) > ch(x) 


c» dist(x, Au ,) 1 


dist(x, r„ 


n + 1 , n + 1 
p ' q 


for an arbitrary value q £ (n + l, oo). Clearly, this would not have sufficed for our 
absorption argument. Thus, the key improvement here consists in the observation 
that G 1 is actually a Lipschitz function away from As the estimates, which 
yield the improved bounds in Proposition 3.10, are interior estimates around A w , 
this regularity away from suffices for our purposes. 


Assuming more regularity on the metric, the lower bound from Lemma 4.13 can 
be further improved: 

Lemma 4.16 (Positivity’). Let a lJ : B + —> be a tensor field that 

satisfies the normalization conditions (A1 )-(Af) from Section 2, condition (ii) from 
Proposition f.6 and in addition is C 1,7 regular for some 7 £ (0,1). Let w : B+ —> R 
be a solution of the thin obstacle problem with metric a 13 and assume that it satisfies 
the normalizations (AO) from Section 2 and (i) from Proposition f.6. Then there 
exists g = rj (n) > 0 such that for e £ C'^ien) 

(76) d e w(x) > cdist(x, A w ) dist(x, r^) - ^. 


Proof. We only sketch the proof here. For each x 0 £ By 2 we consider 
u(x) := d e w(x) + \x' — (xo ) 7 | 2 — 2~ 8 h(x ) — 2 nx 2 +1 , 

where h is as in the proof of Lemma 4.13. Next we split u = zii + U 2 , where u\ 
solves 

dia^djUi — K dist(x, r,„) -2 ui = diF 1 + g, 

where F l = —d e a l ^djU J, g = 2 (xj — (x$)j)dia 1 ^ — Anx n+ idja n+1 '\ 
with K = K(n ) sufficiently large. Hence, u 2 solves 
dia^ djU 2 = g — K dist(x, T w )~ 2 u \, 
where g = -2~ 8 d ia ij d 3 h + 2a u - 4na n+1 ’ n+1 . 

We apply Proposition 3.10 to u\. Since £ C 1 ’ 7 and w £ C 1 ’ 1 / 2 , we have 
dist(x, r w )~ 1 / 2 F t £ L°° and g £ L°°. By Remark 3.12, for any s £ (n + l,oo) 

|ui(x)| < Codist(x, A u ,)dist(x,r u; ) _ ^ + 3 - 

(77) , , n+1 j \ 

• (||F l dist(-,r ro ) 2 |Ix,«(b+) + ||gdist(-, T ^) 1 - 2 II ■ 

For U 2 we apply a comparison argument. By the argument from the proof of 
Proposition 3.17 we obtain that U 2 > 0. Thus u = u\ + 1 x 2 > u\. Evaluating at 
xo yields that d e w(xf) > 2~ 8 h(xo) — tti(xo). Since xq is arbitrary in B.f 2 and by 
using the bound from (77) for u\, we infer that 

(78) d e w(x) > 2~ 8 h{x) — u\(x) > 2 ~ 8 h{x) — c* dist(x, A„,). 

Here we have used that s > 2(n + 1) can be chosen to be sufficiently large and 
||-DcW||l°° < c*. Recalling the estimate for h in (75), we obtain the desired estimate 
for d e w in (76). □ 
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Remark 4.17. We emphasize that without an additional splitting step, we cannot 
hope for an analogous result for a W 1,p , p £ (2 (n + 1), oo], metric a lJ . This is due 
to the fact that by Remark 3.12 for u\ we only have 

|iti(x)| < C 0 dist(x, A w ) p dist(x,P^) 2 - 

• (\\F l dist(-, r w ) 2 ||lp(r«+i) + ||gdist(-, r ^) 1 p 2 ||lp/ 2 (R"+ 1 )) 

< c* dist(x,A UJ ) 1 — dist(x, r^) 2 . 

The bound by dist(x,A u ,) p is optimal in general, due to the divergence right 
hand side diF 1 with F l £ L p . Thus, the resulting bound cannot be absorbed into h~ 
in the estimate (78). 


5. Perturbations 

In this section we consider variants of the variable coefficient thin obstacle prob¬ 
lem. Examples are settings in which the obstacle (c.f. Section 5.2) or the underlying 
manifold (c.f. Section 5.3) is not flat. Moreover, it is also possible to deal with in¬ 
terior thin obstacles (c.f. Section 5.4) and inhomogeneities in the equations (c.f. 
Section 5.1). We show that under suitable conditions on the metric and obstacles, 
it is possible to recover the regularity results from the flat setting. Instead of re¬ 
peating all the necessary arguments, we only point out the main difficulties and 
differences with respect to the flat problem which was discussed in [KRS15] and 
Sections 3-4 of the present article. 

While treating the cases of inhomogeneities, non-flat obstacles and boundaries and 
interior obstacles separately in order to stress the essential aspects of the respective 
situation, we emphasize that it is possible to combine these results into a setting 
involving several/all of these features, e.g. a non-flat obstacle and a non-flat hy¬ 
persurface. 

5.1. Inhomogeneities. Exploiting the scaling of the Carleman inequality, it is 
possible to deal with inhomogeneous thin obstacle problems along similar lines as 
in Sections 2-4 in [KRS15] and Sections 3-4 in this paper. Here the main observation 
is that both for the Carleman estimate and the comparison arguments in Sections 
3 and 4 the inhomogeneity is a “lower order” contribution. 

Proposition 5.1 (Inhomogeneous thin obstacle problem). Let 
a : Bi —» r(”+ 1 ) x ("+ 1 ) be a W 1,p (Bi), p £ (n + l,oo], metric satisfying (Al), 
(A2), (Af). Suppose that f £ L q (B+) for some q £ (n + l,oo]. Assume that 
w £ H 1 (B^) is an L 2 normalized solution to the thin obstacle problem 

dia^djW = f in B+, 
w > 0, —a n+lj i djW > 0, w(a" +1,J djw) = 0 in B[. 

Then, the following statements hold: 

(i) The solution w has the following growth estimate: There exists C > 0 such 
that 

I w(x)\<c{ dist 0 c r ™) 2 (l n (dist(x, P ,,,))) 2 if n + 1 < q < 2 (n + l), 

~ ( dist(x,r u; )i(ln(dist(x,r ! „))) 2 if q>2{n + T). 

The constant C > 0 depends on n,p,q, ||/||l9, ||V<W||lp, ||r’||/ ( 2 ( B +)■ 
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(ii) If additionally £ W 1,P (B+) and f £ L P (B)~) with p £ (2 (n + l),oo], 
then w has the optimal Holder regularity: 

w £ C 1 ^ 2 (B+ 2 ). 

(iii) Under the same assumptions on a 13 and f as in (ii) and assuming that 
0 £ r 3 / 2 (i«), there exist a radius p = p(w,f) > 0 , a parameter a £ (0,1] 
and a C 1,a function g such that (potentially after a rotation) 

r 3 / 2 (wi) nB' p = {x= (x",x n ,0)| x n = g(x”) for x" £ B”}. 

Remark 5.2. Since w solves a linear elliptic equation away from the free boundary, 
we can obtain an interior regularity result for w by combining the growth estimate 
(i) with a local gradient estimate. The proof is standard (see e.g. Proposition 4.8 
in [ 1S15JJ and yields 

|Vw;(a;) — \7w(y)\ 

[ C(q)\x — 2 /| 1 ^“hi 2 (|:r-y|) if n + 1 < q < 2(n + 1), q < p, 

< % C(p,q)\x - if n+l <p < q <2(n+l), 

[ C\x-y\?\n 2 (\x-y\) if min{p, q} > 2(n + 1), 

for all x,y £ B+, 2 . We have 

C(q) —> +00 as q \ 0 and C(p, q) —> +00 as q\p. 

Proof. We only point out the main differences with respect to the previously pre¬ 
sented arguments for / = 0. 


Step 1: Modifications in the Carleman estimate and in the proof of Proposition 
4-5 in [KRS15]. As / £ L q for some q £ (n + l,oo], Uraltseva’s C 1,a regularity 
result remains valid [Ura87]. Furthermore, we remark that it is still possible to 
carry out Uraltseva’s change of coordinates (c.f. Proposition 2.1 in [KRS15] and 
[Ura85]), which transforms the variable coefficient problem into a new variable coef¬ 
ficient thin obstacle problem with an off-diagonal property as in (Al) on B[. Here, 
if / is the original inhomogeneity, the inhomogencity of the transformed equation 
is given by f(y) = f(x)\det(DT(x ))\~ 1 \ x=T -i( y y Consequenty, / inherits the L q 
(and W 1,p ) regularity from /. For convenience of notation, in the sequel we sup¬ 
press the tildas. 

The preceding discussion now allows us to prove the Carleman estimate along the 
same lines as in Section 3 in [KRS15]. The only necessary modification is to incor¬ 
porate a further right hand side contribution: 


(80) 


e T ^|a;| _1 (l + ln(|x|) 2 ) _ U 


< 


CHcq 1 (r 2 C(a ij ) \\e T<j> 


L 2 (A+ r ) 


+ T5 


e T?i (l + ln(|x| ) 2 )-5Vf 


xU 1 w;| 


L 2 (A+ r ) + ll e ^N^IL 2 (A+,.) 


)• 


L 2 (A+ r ) 


where 7 := 1 — and 

' p 


C(a v ) =sup 
At r 


a ij (x) - 6 ij 


+ sup |||x| 7 Va* J | 

p<r<r/2 


L n+ 1 (A+) 
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We estimate the additional term coming from the inhomogeneity for parameters 
t e I 1 ’ i+^r/ 2 ]’ where T o = «o + with k 0 > 0: 


o T <l>\ 


; l/| 


L 2 (A+ r ) 


< 


(81) 


< 


< 


I Li(B+ 




"L*-2(A+ r ) 

r T0 + 1 || Jg. 

"LT**(A+ r ) 


I Lv(B+) ( r P )> ^ — 


= 2 - 


ra + 1 


«0- 


Hence, the contribution originating from the inhomogeneity / in the Carleman 
inequality can be regarded as a “subcritical” contribution (in the sense that <5 > 0) 
for the above range of r if kq < 2 — Analogous to the Carleman estimate, 
Corollary 3.1 in [KRS15] remains valid, if the contribution originating from / is 
included on the right hand side: 


r§(l + | InMD-V^))^- 1 N W 2r , (xo)) 


(82) 


+ rmax{ln(r 2 /ri) ?,ln(r 3 /r 2 ) 1 }e' r ^ (ln(r ' 2 ) ) r 2 1 |MIl=(a+ 2r ( ao)) 


< C(e 


T0(ln(ri)) -1 


kill 


1 


0 Tcf>( ln(r 3 )) -1 


^3 lkHLRA+ i2r3 (x o) ) + 


3 r0(ln(|x|)) 


x \f 


L 2 (A+^ 2r3 (x 0 )) 


)• 


Step 2: Modifications in Section 4 in [KRS15]. Estimate (82) combined with the 
boundedness of (81), immediately allows us to infer the analogue of Proposition 4.1 
in [ ] if k x < 2 — ^±1. Moreover, arguing by contradiction, we also obtain 

that if k x > 2 — ; then not only 


1/2 


In 


lim sup ■ 

r—>-0 


/ w 2 

W+ 


ln(r 


> 2 - 


but also 


1/2 


M / 


lim inf ■ 

r —>-0 


4+ 

i r,2r 


ln(r 


> 2 — 


n + 1 

q 


n+ 1 


Similarly, the doubling property (c.f. Proposition 4.3 in [KRS15]), the blow-up 
result (c.f. Proposition 4.4 in [KRS15]) and the lower semi-continuity statement 
(c.f. Proposition 4.2 in [KRS15]) remain valid at points of vanishing order less than 
2 — Also the growth Lemma 4.1 in [ ] remains valid for all xq £ P^ with 


* = 2-2±i: 


(83) 


sup M < Cr min{ ^o- 2 - £ ^>| ln (r)| 2 , 
Bf (x 0 ) 


where C = C{n,p, ||Va l - 5 || iP(s + ) , ||/|| L9(B +))- 

Apart from these results, we also obtain the existence of homogeneous blow-up 
limits at free boundary points with vanishing order less than 2 — In fact, this 
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is again a consequence of the boundedness of the contributions originating from / 
in the Carleman estimate. 

Therefore, by Proposition 4.4 in [KRS15] and Proposition 4.5 in [KRS15], at each 
x £ Tu, with k x < 2 —tiML, there exists an L 2 normalized blow-up sequence centered 
at x , whose limit is a nontrivial homogeneous global solution (to the constant 
coefficient thin obstacle problem) with homogeneity equal to n x . Then by the 
classification of the homogeneous global solutions (Proposition 4.6 in [KRS15]) on 
the one hand, we obtain that there is no free boundary point with n x < 2 — if 
q £ (n + 1, 2(n + 1)). Hence, (84) turns into: 

(84) sup |ui| < Cr 2 ~-Q -1 ln(r)| 2 , 

B+(x 0 ) 


for all xq £ Tu,. On the other hand, k x > | if q £ [2(n + l),oo]. Combining 
this information with the uniform upper growth estimate (Lemma 4.1 in [KRS15]) 
results in 


(85) 


Cr 2 t (lnr) 2 if q £ (n + 1, 2{n + 1)) 

Cri (In r) 2 if q G [2(n + 1), oo] 


re (0,1/2), 


where the constant C > 0 depends on ||Va y || iP , n,p, ||/|| i? . 


Step 3: Modifications in Sections 3 and f. In order to deduce the regularity of 
the free boundary and to obtain the optimal regularity of w , we argue along the 
lines of Sections 3 and 4. Here we interpret the inhomogeneity in the equation for 
the derivative as a divergence right hand side: 

dia i:t did e w = diF 1 in Bi \ A, 
d e w = 0 in A, 

where for all i £ (1,..., n + 1} 

(86) F i := -(d e a ij )d jW + e i f £ L p (Ri). 

Thus, (potentially) after a rescaling which only depends on ) an d 

II/IIlp(Bi)> we m ay assume that the inhomogeneity is small. Therefore, all the 
results in Sections 3 and 4 are valid. □ 

5.2. Non-flat obstacles. In this section we present the main ideas of dealing 
with non-flat obstacles. Due to our almost scaling critical lower bound in Proposi¬ 
tion 3.17, we are able to deal with the non-flat obstacle problem involving metrics 
a b g 1R 1 ’ P (R 1 ) and obstacles ip £ W 2,P (B 1 ) with p £ (2(n + l),oo]. We however 
stress that the analogues of [KRS15] are valid under the even weaker integrability 
assumption p > n + 1. 

Proposition 5.3 (Non-flat obstacles). Let a ?J : B—> Riym 1 ' )x( '” +1 ' ) be a W 1,p 
metric with p £ (n+l,oo] satisfying (Al), (A2), (Af). Suppose that ip £ W 2,P (B' 1 ). 
Let w : B+ —> R be a solution of the thin obstacle problem 

dia^d^w = 0 in Bf, 

(87) 3 1 

w > ip, —a n+1 ’idjW >0, (w — ip)(a n+1 ’^djw) = 0 on B[. 

Then, the following statements hold: 
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(i) The solution w has the following growth estimate: There exists C > 0 such 
that 

\w(x)\<cl r ») 2 “ ^ r «,))) 2 if n + 1 < p < 2(n + 1), 

~ dist(ic, T^)!(ln(dist(x, T^))) 2 if p>2(n + l). 

The constant C > 0 depends on n,p, |MIw 2 'P(b + ) II u 'IIl 2 (b + )• 

(ii) If additionally a IJ £ W 1,P (B\) and ip £ W 2,P (B' 1 ) with p £ (2(n + 1), oo], 
then the solution w has the optimal Holder regularity: 

w £ C^' 2 {B+ 2 ). 

(iii) Under the conditions in (ii) and assuming that 0 £ T3/2(10), there exist 
a radius p > 0, a parameter a £ (0,1] and a C 1,a function g such that 
(potentially after a rotation) 

r 3 / 2 (w) n B' p = {z = (x",x n , 0)| x n = g(x") for x" £ B"}. 

Remark 5.4. We remark that as in Proposition 5.3 the estimate in (i) immediately 
entails a corresponding regularity result (c.f. Remark 5.2). 

We prove this result by reducing to the setting of flat obstacles. 


Proof. Step 1: Recovery of flat boundary conditions. We first carry out Uraltseva’s 
change of coordinates [ lra85]. Then in order to recover flat boundary conditions, 
we introduce a function v : B—> R with v = w — p. This then leads to an 
inhomogeneous thin obstacle problem with a flat obstacle: 


( 88 ) 


dj/jPdjV = f in B+, 
v > 0, — dn+iv > 0, vd n +\v = 0 on B [, 


with / = —diaPdjp = —(dia i: i)djp — dijp. In recovering the results from the 
flat obstacle setting, the main difficulties arise from the error contributions which 
result from the inhomogeneity. In order to derive (i), we hence argue that the as¬ 
sumptions of Section 5.1 are satisfied. In order to obtain (ii) and (iii) we interpret 
the error as a divergence form right hand side and argue as in Sections 3-4. 


Step 2: Bounding the inhomogeneity. Due to our assumptions on / and a lJ and 
by Sobolev/ Morrey embedding, the metric and the inhomogeneity have the right 
integrability. Indeed, in the setting of (i) we have 
• did** £ L p (B+), a*’ £ L°°(B+), 

. d ZJ p£L p (B[), d jl p£L°°(B' 1 ). 

Hence, / £ L P (B)() which allows us to invoke the results from Section 5.1. This 
then yields the growth estimate stated in (i). 


Step 3: Argument for (ii) and (iii). We argue as in step 3 in Section 5.1. In 
order to obtain the results (ii) and (iii), we consider the equation for tangential 
derivatives of v (after carrying out an odd reflection as described in (15)): 

dial 3 dj e v = - di((d e a lJ )djV ) + d e f in B 1 \ A, 

t 89 -* « n • a 

o e v = 0 m A. 

We interpret the right hand side of (89) as a divergence form contribution with 
(90) F l := — {d e a l ^)djV + e l f, i £ {1,..., n + 1}. 
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The regularity of the metric and obstacle, a u £ W 1 ’ p (Bi), tp £ W 2 ’ P (B' 1 ) with 
p £ (2 (n + 1), oo], entail that 

F l £ L P (B i). 

Moreover, by an appropriate scaling argument (which only depends on || Va 1 - 7 || iP ( Bi ) 
and ||‘/ ? llw 2 'P(B'))’ we ma y assume that is small. Hence, all the results 

in Sections 3-4 remain valid. This then yields the desired results of (ii) and (iii). □ 

5.3. Non-flat boundaries. In this section we very briefly comment on the situa¬ 
tion with non-flat boundaries. In this context we have the following result: 

Proposition 5.5 (Non-flat boundaries). Let Cl C M" +1 be a bounded open subset, 
whose boundary contains the W 2,p hypersurface Ai. Let a 1 - 7 : f l —> Riym 1 ^”^ 
be a uniformly elliptic, symmetric tensor field of class W 1,P (CL), p £ (n + l,oo]. 
Assume that wiHUAf-llisa solution of the thin obstacle problem with zero 
obstacle on the hypersurface Ai: 

dia^djW = 0 in Cl, 

w > 0, Uia^djW > 0, w{via 1 ^djw) = 0 on Ai, 

where u : Ai —> R” +1 denotes the outer unit normal field on Ai. Then, the following 
statements hold: 

1 1— n + 1 

(i) (Optimal regularity) If p £ (n + 1, 2(n + 1)], then w £ C x ' oc v (Cl U Ai); if 
p £ (2(n + 1), oo], then w £ C])^J 2 (Cl U Ai). 

(ii) Assuming that 0 £ r 3 / 2 (nO, there exist a radius p > 0 and a parameter 
a £ (0,1 — such that r 3 / 2 (rc) fl B p is an (n — 1)-dimensional C 1,a 
submanifold. 

Proof. The statement can be immediately reduced to the setting of the usual thin 
obstacle problem by carrying out a change of coordinates (c.f. Proposition 2.1 in 
[KRS15] and [Ura85]) that flattens the free boundary. As the free boundary is W 2,p 
regular, this change of coordinates followed by an application of Uraltseva’s change 
of coordinates from Proposition 2.1 in [KRS15], then implies the integrability and 
differentiability properties (AI), (A2), (A3), (A4). This allows to carry out the 
analysis from [KRS15] and Sections 3-4 of the present article. □ 

5.4. Interior thin obstacles. Finally, we comment on the necessary modification 
steps in obtaining analogous results for interior thin obstacles. In this direction we 
have: 

Proposition 5.6 (Interior thin obstacles). Let a 1 - 7 : B\ -+ be a uni¬ 

formly elliptic, symmetric tensor field of class W 1,p with p £ (n + l,oo], satisfying 
(AI), (A2), (A3) and (A)). Assume that w : B\ —» R is a solution of the interior 
thin obstacle problem with zero obstacle on B(: 

dia^djW = 0 in int B+ U int Bf, 
w > 0, [d n +\w\ > 0, w[d n +iw] = 0 on B(. 

Here [9„+iiu] := (d„+iw;) + — (d n +iw)~ and ( d n +iw ) ± denotes the one-sided trace 
of the fluxes on B(. More precisely, for H n a.e. x £ B( we have 

(d n+ iw) ± (x) := lim d n+1 w(y). 

y^x,±y n+1 >0 
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Then, the following statements hold: 

(i) (Optimal regularity) If p £ (n + 1, 2 (n + 1)], then w £ C io ’ c p (B+ U Bf); 

if p £ (2(n + 1), oo], then w £ U Bf). 

(ii) Assuming that 0 6 r 3 / 2 (u;), there exist a radius p > 0, a parameter a £ 
(0,1 — 2^1] and a C' 1, “ function g such that (potentially after a rotation) 

r 3 / 2 (» n B' p = {x = (x",x n ,$)\ Xn = g(x") for x" £ B p }. 

(iii) There exist functions a,b : r 3 / 2 (w) D B p / 2 —»■ K. with a £ C°’ a (T 3 / 2 (ui) fl 
B p / 2 ) and b £ C' 0 ’5+ a (r 3 /2(wi) D B p / 2 ) for some a > 0 (where a can be 
chosen as in (ii)), such that: for each x 3 £ T 3 / 2 {w) C\B p / 4 , 

— For each e £ S'™ -1 x {0} and x £ B p / 4 (xo), 

3 a(xo)(e ■ v Xo ) / {x - x 0 ) ■ u Xo x n+1 

eW[X) 2 {v X0 • A(x 0 )n Xo y/2 W !/2 ^ {Uxo . A( X o)^ 0 )l/2 ’ {a n+l,n+l {xo)) l/2 

< C(n,p) max{e 0 , c»} dist(x, r u; ) 2 +'>'. 

- For x £ B p / 4 (x 0 ), 
d n +iw(x) - b(x 0 ) 

3 tt(xo) _ ( {x Xq) • I^xq X n4 -\ 

"2(a-+i.«+i(To))i/2 Wl/2 \(n X0 ■ A(x 0 )u Xo y/^ (a"+ 

< C{n,p) max{e 0 , c*} dist(x, r ra ) 2 +T . 

— For x £ B p / 4 (xo) and x — Xq £ sp&n{v Xo , e n+ {\, 

W(X) - i(l0)l “ + ' “ aMWSI2 
< C(n,p) max{e 0 , c»} dist(a;, r^,) 2 + 7 . 

Here w 4 / 2 (x) = Re(x n +ix n+ i) 1/2 , Wi/ 2 (x) = - lva(x n +ix n+1 ) 1/2 , w 3 / 2 {x) = 
Ke(x n +ix n +i) 3 / 2 . Moreover, 7 = min{a, ^±1} and A(xq) = (a* J (a;o)). 

Remark 5.7. As the proof of Proposition 5.6 shows, the function b can be identified 
explicitly as b(x 0 ) = ((9„+iw) + (a;o), where x 3 £ r 3 / 2 (u;) fl B p / 2 . 

Remark 5.8. Similarly as in Sections 5.1-5.3, our results generalize to non-flat 
W 2 ’ p obstacles (p, non-flat W 2,p manifolds A4 and L p inhomoneneities f, where 
p > 2 (n + 1). 

Proof. As in the previous sections, we only comment on the main changes. In the 
first step we derive an analogue of the Carleman inequality (7) in [KRS15]. Since 
an essential ingredient in the proof of the Carleman estimate was the vanishing of 
the boundary contributions due to the complementary boundary conditions, slight 
changes are necessary at this point. As w satisfies an elliptic equation in both the 
upper and the lower half planes, we carry out the Carleman conjugation procedure 
in both half-balls separately. On each of these we obtain boundary contributions 
which do not necessarily vanish. However, adding the contributions originating from 
the lower and upper half-balls allows us to exploit the complementary boundary 
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conditions for the interior thin obstacle problem. Hence, we obtain the vanishing 
of all involved boundary contributions. The Carleman estimate then reads: 

rt(l + | InMD-V^C-))^ 1 M L2(Ar2 2r2{xo)) 

+ Tmax{ln(r 2 /ri)- 1 ,ln(r 3 /r 2 ) _1 }e r ^ ln(, ' 2)) r 2 ' 1 IM| L 2 (Ar2 , 2r2 (x 0 )) 


(91) 


< C(e 


T<t>( ln(n)) -1 


W II 


L 2 (A T 1 : 2 r 1 (x 0 )) 


+ e r0(ln(r 3 )) r -l|| UJ || 


L 2 (A. 


r 3 ,2r 3 


(so)) 


’■0(ln(|:c|))| a .|y 


3 (3:0)) 


where A rir3 (x 0 ) denotes the full annulus around the point Xo £ B[. 

In deriving the analogues of the statements of Section 4 in [KRS15], we modify 
the definition of the blow-up sequences slightly: Indeed, we note that if w(x) is a 
solution of the interior thin obstacle problem, then w £ C 1,a (B+ U B j - ) for some 
a > 0 (c.f. page 207 in [Ura87]). Then for 


the function 

(92) 

solves 


b{x 0 ) := lim d n+1 w{y), 

y->x o ,y n+ i>0 


, {x) := w(x) - b(x 0 )x n+ 1 , 


d l a lJ d j v Xo = f in S 1/2 , 

^xo — 0; [^n+lTxo] — ^xo [^n+lTxo] — 0 On B^j^. 

Here f(x) := — (9ia u )(x)6( xo) £ L p (Bi) by the regularity and integrability as¬ 
sumptions on and the C 1,a (B^ 2 ) regularity of w (c.f. [LJra87]). Therefore, for 
any free boundary point xq £ T^ fl B 1 / 2 the modified functions v Xo (x) now satisfy 

(93) Wxjzo) = 0 = |Vu Xo (x 0 )|. 

Here, by approaching x = Xq from the interior of fl w , we observed that [9 TI _|_iui](a:o) = 
0, which yields the second equality in (93). 

Analogously, we consider the following blow-ups around the point xo, which are 
based on v Xo instead of on w: 

, x V X0 (x 0 + rx) 

V xoA X ) : = _^+l„ „ -■ 

r 2 IM|l 2 (B,.(xo)) 

In the case of the interior thin obstacle problem this replaces the blow-up func¬ 
tions defined in Section 2.2. Arguing on the level of v Xo and v X0:r and using the 
observations from Section 5.1, all the results from Section 4 in [KRS15] follow. In 
particular (93) ensures that the discussion of the vanishing order in Proposition 4.6 
in [KRS15] remains valid. 

Further relying on the functions v Xo , we derive the analogues of the results from 
Sections 3-4 for v Xo by arguing along the same lines as for the boundary thin ob¬ 
stacle problem (in this context, we remark that since solutions are, by definition of 
the interior obstacle problem, already defined in the whole ball, it is not necessary 
to carry out reflection arguments). Hence, all the properties from these sections are 
valid on the level of v Xo ■ These can then be directly transferred to the original func¬ 
tion w. Due to the presence of the normalizing factor b(xo)x n +i this changes the 
asymptotic expansion of d n +iw with respect to the one derived in Remark 4.7 by 
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the constant term b(x o). Using the C 0, “ regularity of a and the triangle inequality, 
we infer from the asymptotics of d n +\w that b £ C 0 , 5 +a (r 3 / 2 (ic) fl B p / 2 ). Indeed, 
for xo,xi £ r 3 / 2 ( 111 ) r\B p / 2 , we fix a point x £ B p with dist^Fu,) ~ \xo — x\\. By 
the triangle inequality we have 

\b(x 0 ) - &(xi)| < |G Xo (x) - G Xl (x)| + C(n,p) max{e 0) c„} dist(x, r^)^, 


where for (er 3 / 2 (w)n B p/2 , 


G ? (x) := - 




2 ( 


n+l,n+l 


( 0) 1/2 


m/2 


(x - £) ■ 


%n +1 


> e • A(O^) 1 / 2 ’ (a-+ 1 .™+i(e)) 1 /2 / 

By using the G 0,a regularity of a, the ellipticity and regularity of A — (a*- 7 ), the 
G 1,a regularity of r 3 / 2 (iiO, as well as the relation dist(x,F^,) ~ |x 0 — Xi|, we have 


|G X 0 (x) - G xi (x )| < \x 0 - xi| = 


+ Q! 


This proves the G 0 , 2 + “ regularity of b. The remaining observations follow along 
similar lines as in Sections 3-4. □ 
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